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As  part  of  an  investigation  of  the  physical  basis  underlying  aspects  of  MOS 
processing  technology,  we  have  studied  the  transport  mechanism  by  which  water 
diffuses  through  thin  Si02  films.  This  process  is  responsible  for  the  formation 
of  oxide  layers  on  silicon  wafers  by  means  of  thermal  steam  oxidation,  which  is 
frequently  employed  in  the  manufacture  of  integrated  circuit  devices.  We  have 
performed  tracer  diffusion  measurements  involving  network  , demonstrating  the 

importance  of  oxygen  exchange  between  the  Si02  network  and  molecularly  dissolved 
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water „  We  have  foun^that  in  the  presence  of  water,  bound  network  oxygen 
diffuses  through  SiC>2  as  a  constituent  of  molecularly  dissolved  water.  Employing 
methods  common  to  state-of-the-art  semiconductor  technology,  the  central  region 
within  a  thermal  oxide  layer  grown  on  silicon  was  enriched  with  immobile  by 
ion  implantation.  After  heating  in  atmospheres  with  different  water  contents, 
the  extent  of  diffusion  was  determined  by  observing  changes  in  the  concen¬ 
tration  profile  (i.e.  the  chemical  concentration  as  a  function  of  depth)  of 
implanted  by  means  of  nuclear  resonance  profiling,  .  In  this  technique,  high- 
energy  protons  from  a  Van  De  Graaff  accelerator  lose  energy  in  their  passage 
through  the  oxide  until  their  energy  falls  to  a  resonant  energy  of  the 
i  18o(p,«)^^N  reaction.  High-energy  alpha  particles,  emitted  in  this  inelastic 
I  reaction  from  the  which  is  located  at  that  depth,  are  detected  after 
'  escaping  from  the  surface.  The  observed  thermal  broadening  of  ion-implant 

b ■ stributions  permitted  bulk  diffusivity  to  be  measured  conveniently  at  temper- 
•  •ures  or  gas-phase  concentrations  much  lower  than  previously  possible.  Dif- 
■s ions  conducted  in  steam  at  1  atm  (at  temperatures  as  low  as  250C)  showed  an 
r  t I vat ion  energy  of  about  17  kcal/mol,  which  is  close  to  that  originally 
•’'ensured  for  water  permeation  in  Si02.  Diffusions  in  both  air  and  dry  nitrogen 
showed  a  similar  activation  energy,  with  respective  ore-exponential  factors  two 
■  J  three  orders  of  magnitude  below  the  steam  value.  Diffusions  in  low  pressure 
v  'or  vapor  showe  '  a  clearly  linear  dependence  on  gas  phase  water  concentration 
•vu  to  80  ppm,  in  contrast  to  thermal  oxidation  results  reported  by  others 
The  results  of  this  study  are  consistent  with  a  recently  proposed  model  of 
w'  ter  diffusion  in  Si02, in  which  the  diffusion  mechanism  is  the  interstitial 
-Tvinsport  of  dissolved  molecular  water  accompanied  by  a  reversible  reaction  with 
: icon-oxygen  bonds  in  the  network. 
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NUCLEAR  TRACER  MEASUREMENTS 
OF  LOW  TEMPERATURE  WATER  DIFFUSION 
IN  Si00  THIN  FILMS 


INTRODUCTION 


Background 

Ttie  formation  of  thin  insulating  films  on  the  surface  of  semiconductor 
materials  is  an  important  part  of  modern  integrated  circuit  technology,  which 
has  assumed  a  central  role  in  the  production  of  Army  control,  communications, 
radar  and  data  processing  systems.  The  thin  film  can  serve  in  integrated 
circuits  as  a  dopant  mask,  surface  protector,  primary  passivation  coating, 
component  isoLator,  and  as  gate  insulacor  in  MOS  structures  (1).  Other 
applications  of  current  interest  include  solar  cells  and  power  transistors. 
Because  of  their  military  and  commercial  importance,  various  means  of  film 
growth  have  been  intensively  studied,  especially  for  high-quality  gate 
--ins.ulators_on  silicon  surfaces.  For  that  application  a  more  or  less  standard 
practice  of  thermal  oxfdalTon  1  ri~ cTi*y~ y  g  ott  o  t"  s  t e-n^c-a-c-ms-s-ph a  r-s-x- -hi -s-  .s-UQ.Lveri.-j. 
A  properly  made  thermal  oxide  film  must  have  a  high  dielectric  strength  and 
immunity  to  electrical  breakdown,  contain  an  acceptably  low  amount  of  trapped 
charge  (particularly  at  the  oxide-silicon  interface),  be  uniform  in  thickness, 
bond  intimutelv  to  the  semiconductor  substrate,  and  resist  chemical  or  mechan¬ 
ical  degradation.  These  requirements  are  well  satisfied  by  present  thermal 
oxidation  methods  However,  as  integrated  circuit  devices  become  more 
complex  ana  operate  at  higher  speeds,  they  require  smaller  and  more  closely 
packed  component  structures.  The  scale  of  sizes  being  contemplated  has  now 
reached  a  point  where  basic  changes  will  soon  be  required  in  manufacturing 
methods.  For  example,  the  need  is  foreseen  to  hold  temperatures  in  all  steps 
of  a  processing  sequence  low  enough  to  prevent  significant  migration  of 
dopants.  This  has  spurred  work  on  low- temperature,  high-pressure  thermal 
oxidation  (2,3). 

As  might  be  expected,  these  low-temperature  processing  techniques  may 
offer  enough  advantages  over  present  high-temperature  processing  to  supplant 
it  ir.  the  manufacturing  technology  of  future  integrated  circuit  devices. 

Tiio.se  techniques  may  well  produce  devices  having  electrical  characteristics, 
lifetimes  and  reliabilities  surpassing  any  previously  attainable.  The  ubi¬ 
quity  of  Si02  films  oh  semiconductor  materials  has  been  accompanied  by  a 
vast  research  effort  -  a  recent  bibliography  (4)  lists  560  references. 

Included  are  many  studies  of  film  formation,  particularly  thermal  oxidation. 


1.  A.  S.  drove,  Phys ics  and  Technol ogy  o f  Semiconductor  Devices (Wiley,  1967) 

2.  K.  i.  Zeto,  C.  Thornton,  E.  Hrvckowian  and  C.  D.  Bosco,J  Kleetrochem 
Soc  112,  1411  (1975) 

3.  L.  E.  K.iz  anu  i».  F.  Howells,  .[  Electrochem  Soc  126,  1822  (1979) 

4.  H .  Ago  j.xnian,  Sol  St  Tech  19,  Jh  (1977) 
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V  .  l  tualiy  all  of  this  work  has  been  device-oriented.  Physical  models  of  the 
oxidation  process  have  received  little  attention;  only  now  are  models  being 
tested  with  sufficient  care  to  generate  a  consensus  of  acceptance.  Through¬ 
out  the  device-oriented  literature,  one  encounters  a  virtually  unanimous 
lack  of  appreciation  of  both  the  nature  and  the  consequences  of  the  water- 
Si02  reaction.  The  lack  of  reproducibility  in  thermal  growth  kinetics  data 
caused  by  the  presence  of  water,  even  in  trace  amounts,  has  lately  been 
corrected  through  careful  control  of  conditions  (5).  For  the  future,  a 
greater  understanding  of  the  mechanism  responsible  for  oxidation  in  the 
presence  of  water  is  imperative,  especially  in  the  new  low-temperature  regime 
where  it  now  appears  that  processing  of  the  coming  generation  of  devices  will 
.  tk  place. 


"he  purpose  of  this  investigation  was  threefold: 

"o  clarify  the  mechanism  of  SiC>2  film  growth  in  atmospheres  contain- 

_ through  its  determination  of  the  kinetics  of  network  oxygen 

i  ,  this  work  provides  the  First  3 Free tTT’xpe r imentaT  demons tfr a  t  Tor  oF 
•  ot  a  recently  proposed  model  of  the  oxidation  process.,  Aceord- 

■  this  model,  growth  proceeds  by  inward  interstitial  diffusion  of 
'ved  molecular  water,  which  simultaneous i y  reacts  with  the  SiOy  network 
immobile  reaction  products. 

!’o  apply  the  technique  of  ion  implantation  to  measure  tracer  oxygen 
>n  in  Si02.  The  observation  of  tracer  distribution  broadening  both 
.as  surface  effects  and  enhances  the  sensitivity  of  diffusion  measure- 
that  tracer  diffusion  measurements  can  be  carried  out  at  eonsidera- 
.'red  temperatures. 

To  provide  a  framework  for  resolving  several  long-standing  discrepan- 
’ ".'tween  measurements  involving  water  and  oxygen  diffusion  in  Si02. 

"■;,v  .ire  wide  d  i  f  f  ereneos  in 'oxygen  di  l  I'usivi  t  ies  observed  in  perinea - 
•  ■-.ner  imenLs,  as  well  as  in  measurements  involving  tlie  pressure  depend- 
of  the  growth  kinetics  of  thermal  silicon  oxidation^  This  work  demon- 
that  these  can  be  resolved  through  a  proper  understanding  of  water 
■  :-.Lor>  in  tiie  oxides. 


ill  Structure  of  S10 ? 

basic  building  block  of  vitreous  silica,  the  amorphous  form  of  SiC^, 
'■t rahed ron  of  oxygen  ions  around  a  central  silicon  ion.  According  to 
•’•aPy  accepted  model  originated  by  Zaehariasen  (6),  each  tetrahedron 
■■ted  to  four  others  by  silicon-oxygen  bonds  at  its  corners  in  a  con- 
• hreo  dimensional  random  network  lacking  any  symmetry  or  periodicity. 


.  S  Mayo  and  W.H.  Evans,  "Development  of  Hydrogen  and  Hydroxyl  Contamina¬ 
tion  in  Thin  Silicon  Dioxide  Films"  NBSIR  78—1558  (NBS,  March  1979) 

ft.  W.  H.  Zaehariasen,  .1  Am  Chem  Sob  54  3841  (1932) 
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as  illustrated  in  Figure  i .  Thu  Si-O  interatomic  distance  is  1.62  A;  the 
0-0  distance  is  2. 65  A.  The  distribution  of  Si-O-Si  bond  angles  ranges  from 
120°  to  180°;  it  is  skewed  toward  the  smaller  angles  and  has  a  single  maximum 
located  at  144°  (7).  It  is  narrow  compared  to  a  fully  random  distribution, 
so  that  the  structure  is  quite  uniform  over  several  tetrahedra  (although 
not  beyond  that).  There  is  a  good  deal  of  experimental  evidence  supporting 
this  model  which  is  quoted  by  Doremus,  pp.  24-28  (8),  including  X-ray  diffrac¬ 
tion  studies  in  terms  of  pair  distribet ion  functions,  transmission  electron 
microscopy  and  atomic  vibrational  spectroscopy  via  infrared  absorption,  Hainan 
emission  and  inelastic  neutron  scattering.  Beside  this  work  are  studies 
involving  X-ray  diffraction  (9,  10),  inelastic  neutron  scattering  (11),  and 
electron  diffraction  and  infrared  absorption  (12). 

Several  recent  experiments  (13,  14,  15)  support  closely  related  but  com¬ 
peting  models  based  on  a  collection  of  distorted  microcrystallites  such  as 
crystpbalite,  each  containing  only  a  few  tetrahedra.  All  viable  models, 
however,  share  the  essential  picture  of  vitreous  silica  as  a  loosely  packed 
material  in  which  cavities  or  interstices  are  surrounded  by  walls  of  more  or 
less  intact  SiOj.  tetrahedra.  Shackelford  and  Masaryk  (15)  give  the  mean 
diameter  of  the  interstices  as  1.96  X. 

There  is  much  evidence  that  ary  thermal  oxide  layers  formed  on  silicon 
surfaces  have  Lae  same  structure  as  vitreous  silica  X-ray  diffraction 
experiments  showed  tee  oxide  layers  to  oe  amorphous  (16).  The  density  of 
the  lay>.r-  is  c  .use  to  tnat  of  fused  silica  (17),  as  are  the  changes  in 
density  emteti  by  ionizing  radiation  (18).  The  optical  properties  are 
essentially  identical  (19),  as  are  the  Si-0-Si  bond  angles  obtained  from 
electron  diffraction  and  JR  absorption  experiments  (12). 


7.  X.  ;  .  Mo/.z  and  B.  E.  Warren,  .1  A. .pi  Cryst  2,  164  (1969) 

3.  R.  H.  i)orem.,.s ,  GJ ,i_ss  5.  ier.ee  (Wiley,  1973) 

9.  J.  H.  Konnert  and  a.  Karie,  Acta  Cryst  A  29,  702  (1973) 

10.  A.  :i ,  Marten,  J  Chum  Fays  5_f>,  1905  (  1972) 

11.  A.  L.  headbutt or  and  M.  W.  Stringfellow,  "Neutron  Inelastic  Scattering," 
Proc  Grenoble  Conf  (IAEA,  Vienna,  1972)  p.  501 

12.  X.  Nagas ima ,  Japan  J  Appl  Phys  9,  879  (1970) 

13.  0.  Bratter,  C.  fioden,  A.  Balogh  and  A.  Andreel:,  Appl  I'hvs  16,  211 
( 19/8) 


.14.  P,  H.  Gaskell,  1).  W.  Johnson,  J.  Mumrvst  Sol  20,  1/1  (19/6) 

15.  J.  !•'.  Shackelford  and  J.  S.  Masai  vk,  ’  Moneryst  Sol  30,  127  (1978) 

16.  M.  B.  Brodsky,  I).  Cubicciotti,  J  Am  Chem  Soc  7_3,  3497  (1951) 

17.  R.  J,  Doremus,  J  Phys  Chem  80,  1773  (1976) 

18.  K.  P.  KerNisse,  J  Appi  Phys  4_5,  167  (1974) 

19.  A.  N.  Knopp  and  R.  Stickler,  Electrochem  Technology  _5 ,  37  (1967) 
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Ion  tmplaiiLat  ion  .nnl  i  i s  K fleets  on  s  i  q  , 

The  implantation  of  heavy  ions  in  the  tens  to  hundreds  of  keV  range  is 
by  now  a  routine  aspect  of  semiconductor  processing  (20).  The  implantation 
of  dopant  elements  directly  into  desired  regions  in  the  interior  of  semi¬ 
conductor  materials  has  eliminated  the  need  for  high-temperature  diffusions 
to  induce  their  migration  from  the  surface,  allowing  much  greater  control 
over  material  properties.  In  regular  lattice  structures,  the  ions  produce 
heavy  damage,  and  a  significant  fraction  of  them  remain  interstitial.  When 
appropriate  annealing  procedures  are  employed,  the  damage  is  healed  virtually 
completely  and  the  dopant  ions  assume  desired  substitutional  sites  (21). 

There  have  been  many  studies  of  ionizing  radiation  effects  on  SiO^  as 
well.  One  motivation  for  these  is  nuclear  vulnerability  assessment  for 
military  applications;  another  comes  from  tne  frequent  practice  of  implan¬ 
tation  into  silicon  through  previously  grown  Si02  layers.  Many  different 
effects  have  been  observed,  including  radiation  compaction  (22,  23,  24,  25, 
18),  enhancement  of  HF  etch  rate  (26,  27),  accumulation  of  thermal  energy 
(28),  changes  in  Raman  and  IR  spectra  (29,  30),  change  in  refractive  index 


20.  J.  W.  Mayer,  L.  Eriksson  and  J.  A.  Davies,  Ion  Implantation  in  Semi¬ 
conductors  (Academic,  1970) 

21.  J.  R.  Gibbons,  Proc  IEEE  60,  1062  (1972) 

22.  W.  Pr.unak ,  J  Appl  Phys  4_3,  2745  (1972) 

23.  J.  E.  Shelby,  J  Appl  Phys  _50,  3702  (1979) 

2q.  C.  8  Norris  and  E.  P.  EerNisse,  J  Appl  Phys  45,  3876  (1974) 

25.  E.  P.  EerNisse  and  C.  8.  Norris,  J  Appl  Phys  U5_,  5196  (1974) 

26  W.  Rratschmer,  "Effects  of  Heavy  Ion  Radiation  on  Quartz  Glass,1' 

Prni  .  lot.  l.unf.  m)  Min’  I  Photography  and  Track  Detectors 
(Bucharest  L9/2)  quoted  in  Antonini  (1978) 

27.  A.  Nonfret  and  J.  Bernard,  "Chemical  and  Electrical  Behavior  of  Ion 
Implanted  S107  Films”  in  Proc.  2nd  Int.  Conf.  Ion  Implantation  in 
Semiconductors,  I.  Ruge  and  J.  Graul,  eds.  (Springer-Verlag,  1971). 

28.  V.  Antonini,  A.  Manara  and  P.  Lensi,  "Ion  Irradiation  and  Stored  Energy 
in  Vitreous  Silica,  in  Pantelides  (1978),  p.  316 

29.  J.  R.  Bates,  R  W.  Hendricks  and  L.  B.  Shaffer,  J  Chem  Phys  6^,  4163 
(1974) 

30.  C.  K.  Frxtzcne  anu  W„  Rothermund,  J  Electrochem  Soc  119 ,  1243  (1972) 
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(ibid),  thermo  luminescence  (31),  and  generation  of  deep  electron  traps 
33  ,  34)  and  color  centers  (35,  36,  37). 


The  mechanisms  responsible  for  these  effects  fall  into  three  rough 
categories,  which  are  separable  to  some  extent  by  choice  of  appropriate 
projectile  and  by  studying  the  annealing  behavior  of  the  defects  created. 

The  first  is  bond  breaking  resulting  from  ionization  reactions,  exemplified 
by  dangling  bonds  and  oxygen  vacancies.  These  are  produced  by  X-  and  gamma 
radiation  and  electrons,  and  are  largely  annealed  out  by  heating  to  400  C; 
some  remain  until  700  C  (25) .  The  second  is  formation  of  new  bonds  between 
implanted  heavy  ions  and  the  network.  These  do  not  anneal  out  but  remain  as 
features  of  an  altered  network  (34).  The  third,  structural  alterations 
caused  by  atomic  displacement,  is  associated  with  heavy  ion  bombardment 
These  may  be  produced  by  direct  nuclear  knock-on  reactions  causing  displace- 
■fill  ot  recoils,  as  we  I  I  as  through  absorption  of  sufficient  vibrational 
energy  to  reconfigure  strained  bond  structures.  Their  importance  increases 
with  increasing  ion  mass  and  decreasing  velocity,  being  especially  signifi- 
■  ■  -i t  near  the  end  of  the  ion  paths;  most  displacement  damage  is  confined 
■bat  region  (38,  39).  An  additional  effect  peculiar  to  the  presence  of 
'•••'non  has  recently  been  identified  by  Shelby  (23)  as  dilatation  or 

leg  of  the  network,  possibly  caused  by  a  radiation- induced  increase 
•.be  number  of  binding  sites  for  SiOH  and  SiH.  Irene  and  Chez  (401  have 
-euJated  that  this  dilatation  may  contribute  to  the  pronounced  effect  of 
•  o  .  II  amounts  of  water  in  the  network. 

As  in  the  case  of  silicon,  inert-atmosphere  heating  has  been  observed  to 
>••  'en.'  out  the  various  radiation- induced  alterations  to  the  Si02  network. 


T.  Shelby,  J  Appl  Phys  50,  3702  (1979) 

'•  P.  EerSisse  and  C.  B.  Norris,  .1  Appl  Phys  4_5,  5196  (1974) 

3’.  G.  W.  Arnold,  "Thermoluminescence  in  Ion-Implanted  Si07,"  in  Chernow 
et  al  (1977) 

.  P.  Donovan  and  M.  Simons,  J  Appl  Phys  43,  2897  (1972) 

Johnson,  W.  C.  Johnson  and  M.  A.  Lampert,  J  Appl  Phvs  4_6,  1216 

C  9  7  I) 

'*•  •  l)i  Maria,  1).  R.  Young,  W.  R.  Hunter  and  C.  M.  Serrano,  IBM  J  Res 

bevel ,  22,  289  (1978) 

Arnold,  IEEE  Trans  Nucl  Sci  NS-20,  220  (1973) 

.  Arnold,  "Vibrational  and  Electronic  Spectroscopy  of  Ton-lmplanta- 
■  on- !  ndeeed  Detects  in  Fused  Silica  and  Crystalline  Quartz,'.'  in 
ic.-s,  p.  ,'/K  M978) 

P •  Sige! ,  Jr.,  B.  I).  Evans,  R.  J.  Ginther,  E.  J.  Friebole,  D.  L. 

'■rt  scorn  and  NRI.  Memo  Rpt  No.  2934  (NKL,  1974) 

38.  D.  K.  Brice,  Radiation  Eff  (OB)  6,  77  (1970) 

39.  K.  B.  Winterbon,  Ion  Implantation  Range  and  Energy  Depositions,  Vol  2 
(Plenum,  1975) 


40.  K.  A.  Irene  and  R.  Chez,  J  Electrochem  Roc  124,  1757  (1977) 


Antonin i  t*e  ai  (.2d)  observed  chat  the  t  hcrmai  energy  ..ctuiwJtiot  during 
bombardment  with  1  MeV  alpha  particles  was  released  upon  heating  to  about 
470  C,  while  that  accumulated  during,  bomba rdme tt  with  46.5  MeV  Ni+6  ions  was 
releasee  in  two  stages:  partly  at  450  C  anu  tne  rest  at  about  610  C. 

Monfret  and  Bernard  (27)  reported  that  radiation  enhancement  of  the  HF  etch 
rate  largely  disappeared  after  a  20  min  anneal  at  600  C  (totally  after 
1000  C) .  Wang  et  al  (41)  reported  tne  changes  (induced  by  50  keV  A1  ion 
implantation)  in  i’V  absorption  in  the  7.5  eV  band  to  be  7  5%  gone  with  a 
500  C  anneal  (totally  after  900  C) ,  and  that  the  amplitude  of  a  similarly 
i:\duc oe  HSR  resonance  (with  g  =  2.0037)  decreased  with  increasing  tempera¬ 
ture.,  disappearing  by  500  C.  Eernisse  and  Norris  (25)  in  their  thorough 
siuey  of  compaction  annealing,  reported  a  relatively  sharp  return  to  initial 
density  as  temperatures  increased  beyond  650  C  after  bombardment  by  either 
7.27  keV  electrons  or  250  keV  protons,  and  a  gradual  return  between  300  C 
and  900  C  after  500  keV^®  Ar'r+  bombardment.  In  sum,  and  as  indicated  by 
Sige.l  (42),  all  studies  of  radiation  carnage  in  SiOy  indicate  that  a  900  C 
anneal  is  sufficient  to  completely  restore  the  network  to  its  pre-irradiation 
s  ta  Le «. 

The  be.iavior  of  the  implanted  ions  at  the  end  of  their  range  in  the  net¬ 
work  car.  be  mi  -.-i  r  to  using  a  model  •  f  ion-network  interactions  propounded  by 
Frimak  197  2  (22/  ana  Pritaak  1975,  p>.  7.11-119  (4  3).  This  is  a  refinement  of 
the  c hermit!  spike  model  trig: naliy  proposed  by  Seitz  (44),  in  which  heavy 
ions  transfer  energy  to  the  network  vs...  atomic,  collisions,  causing  a  large 
amoj'it  of  vibrational  excitation.  Subsequent,  rapid  refreezing  leads  to 
reformat ion  of  fco.vus  with  the  implanted  ions  in  substitional  sites.  Evi¬ 
ct  once  ‘ram  this  comes  from  SIMS  measurements  by  H.  Hughes  (quoted  in  Johnson 
et  al  (33)),  who  found  that  30-min.  anneals  in  N'2  at  900  C  failed  to  alter 
the  concentration  profiles  of  Al  implanted  in  Si02;  a  similar  result  was 
obtained  by  DiMaria  et  al  (34)  at  1050  C. 

22.  W.  Prinak,  J  App.L  Phys  43,  2745  (1972) 

25.  K.  P.  HerN.sse  ana  C.  6.  Norris,  J  Appl  Phys  45,  5196  (1974) 

27,  A.  Mc.nfrec  and  J.  Bernard,  "Chemical  and  Electrical  Behavior  of  Ion 
Implanted  Si02  Films"  in  Proc,  2nd  Int .  Conf.  Ion  Implantation  in 
Semiconductors ,  I.  Ruge  and  J.  Graul,  eds.  (Springer-Verlag,  1971) 

2b.  V.  Antonin! ,  A.  Manara  and  P.  Lensi,  "Ion  Irradiation  and  Stored  Energy 
in  Vitreous  Silica," in  Pantelides,  p  316  (1978) 

33.  X.  M  Johnson,  W.  C.  Johnson  and  M.  A,  Lampert,  J  Appl  Phys  4_6,  1216 

(1975) 

34.  D.  J.  biMaria,  D.  R.  Young,  W.  R.  Hunter  and  C.  M.  Serrano,  IBM  J  Res 
bevel,  22,  289  (1978) 

41.  S.  Wang,  T.  Russell  and  B.  S.  H.  Royce  ^'Annealing  Studies  of  A1+ 
Implanted  Si02  Thin  Films."  PSSL  300874  (Princeton  U.,  1974) 

42.  G.  H.  Sigel,  private  communication  (1.979) 

4  3 .  W„  ,’ri.mak.  Compacted  States  of  Vitreous  Silica  (Gordon  and  Breach,  1975) 
44.  .'.  Seitz,  Disc  Faraday  Soc.  5_,  271  (1949) 


7 


' r a n  sp < i  r  t  u f  0 and  H,,0  in  SiO^ 


Definitions  of  Solubility 

The  solubility  generally  found  in  the  literature  is  the  saturation 
concentration  of  gas  dissolved  in  the  solidyC>dissolved (sa t )  (in  units  ot 
molecuies/cm^)  at  standard  temperature  and  pressure.  A  related  quantity, 

•  :a  solubility  ratio,  also  frequently  appears  in  the  literature.  This 
quantity  is  defined  as  the  steady  state  ratio_^jissoiveci(sat)  //?^as  phase 
;  ,i  given  temperature  and  pressure.  This  is  a  dimensionless  quantity.  It 
,,  U  pi  iHlenl  ot  ambient  gas— phase  pressure  provided  the  same  substance  is 
in  observed  both  outside  and  inside  the  oxide  (l.ts  in  the  absence  o( 
d  ‘  h:~oc  i  at  ion) .  Its  temperature  dependence  -  generally  very  weak  -  arises 

possible  temperature  dependence  in  the  activation  energy  oi  molecular 
I  4  5) 


discussed  shortly,  the  dissolution  of  wate 


sm 


?  involved 


reaction  with  the  network,  in  which  the  molecuiarly  dissolvec  water 
Lutes  to  form  two  OH  units  with  the  additional  0  ion  being  provided  by 
network  (4b,  45).  This  reversible  reaction 


H20  +  Si-O-Si  2  Si OH 


(1) 


e  equilibrium  concentrations  of  reaction  products  which  are  large 
to  that  of  molecuiarly  dissolved  water.  The  equilibrium  concent  ra¬ 
re  related  by  the  law  of  detailed  balance 


K 


2  >^0 H 

W  * 


(2) 


2 

o  or  K  (including  any  temperature  dependence)  is  gotten  from  the 
’  ->1.11  ties  of  molecular  and  "reacted"  water  (i.e.  SiOH  groups).  It 

•  ■.■member ed  that  the  dissociation  of  H20  leads  to  a  reacted  water 
.v  which  is  dependent  on  the  concentration  of  vapor  phase  water. 


4).,  .! .  F.  Shackelford  and  J.  S.  Masaryk,  J  Noncryst  Sol  2J_,  55  (1976) 

46.  R.  H.  Doremus,  in  Reactivity  of  Solids,  Mitchell,  de  Vries, 
Roberts  and  Cannon,  Eds.,  0-ri lev,  In6°)  p.  667 
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Definitions  <>;  Diffusivitv 


fhere  are  three  difierent  diffusion  roe:  i  icients  which  are  involved  in 
the  transport  of  tracer  The  first  of  these,  Jeff,  is  that  quantity 

normally  measured  in  permeation  experiments.  The  basic  method  of  determining 
this  effective  diffusivity  is  by  observing  the  steady-state  flux  J  of  gas 
passing  through  a  membrane  of  thickness  X0  whose  two  surfaces  are  main¬ 
tained  at  different  pressures.  From  the  basic  equation  J  =  -Dt,f 
we  have 


issolved 


el  r 


_Xo _ 

^Pd  issolved 


(3) 


where  0_  ■>.  , ,  the  aifference  in  concentrations  at  the  two  surfaces,  is 

known  through  independent  measurement  of  the  solubility.  The  temperature 
dependence  of  1\.  r\  is  generally  expressed  via  the  Arrhenius  equation 


Ea  \ 


(4) 


where  K  and  tC  ar-.-  the  universal  gas  constant  (in  units  of  kcal/mol/  K  and 
eV/  k  res. .actively)  uc.u  0  Ea  are  tne  activation  energy  of  diffusion  (in 
units  ci  kca...nui  ana  eV  respectively).  Jeff  can  also  be  measured  dynami¬ 
cally  ay  et serving  the  evolution  of  ambient  gas  pressures  in  the  early  stages 
of  eitner  permeation  or  desorption  in  membranes,  filaments  or  powders  in  an 
enclosed,  volume  (see,  for  example.  Dor emus,  pp.  128-130  (8)),  It  should  be 
no tee  caul  Jeff  will  be  independent  of  ambient  gas-phase  concentration  only 
if  till  fluxes  of  both  the  emerging  particles  and  the  dissolved  species  have 
tne  same  pressure  dependence.  Such  is  not  the  case  for  water  diffusion, 
where  De- ;  is  the  apparent  diffusivity  of  reacted  water. 


£k  transport  of  the  actual  nooiie  species  is  described  by  the  micro¬ 
scopic  diffusivity  0.  It  is  influenced  only  by  the  interactions  of  individ¬ 
ual  diffusing  molecules  with  tne  network.  For  a  dilute  gas  (i.e.  one  with¬ 
out  significant  interaction  between  diffusing  molecules)  in  a  spatially 
uniform  network,  D  is  independent  of  ambient  gas  phase  concentration, 
spatial  Location  and  dissolved  gas  concentration  -  it  depends  only  on  temper¬ 
ature.  For  oxygen  diffusion  in  SiC>2,  which  proceeds  without  any  appreciable 
dissociation  or  reaction  with  the  network,  Deff  and  D  have  the  same  value. 
However,  for  water  diffusion,  they  are  related  by  the  equation 


Deff  =  D 


2/?H2 


2y0H2O  + p 


OH 


(5) 


the  fraction  represents  the  fraction  of  rime  that  any  particular  OH  group 
finds  itself  moving  as  part  of  a  diffusing  water  molecule.  The  factor  of  2 
appears  because  each  water  molecule  removes  two  OH  groups  from  the  network. 


i; .  II.  a.irenms,  OI.i.ss  Science  (Wilev,  1071) 
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i  he  third  diffusion  coefficient,  D*,  describes  the  apparent  diffusion  of 
tracer  ^0.  As  wiLl  be  discussed  shortly,  the  mechanism  responsible  for  the 
m  Lizr.it  ion  of  network  -^0  observed  in  this  study  is  the  exchange  between 
difi using  water  and  network  oxygen.  The  reaction  between  dissolved  molecular 
water  and  the  SiC>2  network  described  by  equation  (1)  acts  to  remove  from 
its  network  site  and  to  convert  it  instead  to  a  constituent  of  diffusing 
molecular  water: 

H2160  +  Si-l80-Si  ->  Si160H  +  Si180H  ->  H^O  +  Si-160-Si  . 


•'  the  presence  of  free  exchange,  the  tracer  diffusivity  D*  is  related  to  tin. 
< i t usivity  of  molecularly  dissolved  water  by  the  equation 


D*  =  D 


A  2° 

+/?Ho0 


lion  is  the  fraction  of  time  that  an  J  0  finds  itself  being  carried 
■  a  diffusing  water  molecule.  The  factor  of  2  accounts  for  the  two 
.  oms  per  Si02  molecule. 


Mechanisms  of  02  and  H20  Diffusion 

vised  on  the  observations  to  be  discussed  shortly,  the  mechanism  of 
•  oxygen  or  water  diffusion  has  been  taken  to  be  the  interstitial  trans- 
-  dissolved  molecular  oxygen  or  water  respectively.  In  neither  case 
■’-e  diffusion  pro. ess  involve  direct  bonding  between  the  diffusing 

and  the  ,S102  network.  Rather  it  involves  the  transport  of  nonreacting 
■■.'ey,  through  the  interstices  or  cavities  in  the  network.  This  inter¬ 
model  was  first  proposed  by  Anderson  and  Stuart  (47);  they  calculated 

■  ! vat  ion  energy  for  diffusion  as  the  elastic  energy  required  to  deform 

■  by  en’argrng  the  circular  doorwavs  between  interstices  enough  for 
’.'using  species  to  pass.  An  essential  feature  of  this  model  is  a  lack 

•icmicaJ  binding  between  the  diffusing  species  and  the  network  (Milberg 
1  S' ) .  Although  the  model  is  simplistic,  the  functional  dependence  of  activn- 
■nergy  upon  molecular  diameter  agrees  well  with  activation  energies 

for  molecular  diffusion  of  various  gases  in  fused  silica  (Doremus, 
r S ) > .  It  is  also  consistent  with  the  linear  dependence  of  diffusivity 
■"’•lont  pressure  observed  bv  Norton  (49)  for  oxygen,  whose  diffusion  is 
umpanied  by  any  significant  reactions  with  the  Si02  network. 

Doremus,  Class  Science  (Wiley,  1973) 

('  h.  Anderson  and  D.  A.  Stuart,  J  Am  Ceram  Soc  3J,  573  (1964) 

n3.  M.  E.  Milberg,  "Diffusion  in  Glass,"  in  Fast  Ion  Transport  In  Solids, 

W.  Van  Gool  ed.  (North  Holland  1973)  p.  37S 

49.  F.  .1.  Norton,  Natuie  19_1,  701  (1961) 
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inis  was  demons tratea  recently  in  a  tracer  experiment  conducted  by  Rosencher 
et  al  (50)  in  which  thermal  oxide  films  which  were  first  grown  in  natural 
oxygen  were  further  grown  in  loO-enriched  gas.  Almost  all  additional 
was  observed  to  collect  in  a  layer  at  the  si— Si02  interface,  while  less  than 
0.1%  remained  in  the  previously  grown  region.  Great  care  was  taken  to 
eliminate  any  water  in  this  experiment. 

The  mechanism  for  water  diffusion  tnrough  SiOp  is  similarly  thought  to 
involve  the  interstitial  transport  of  dissolved  molecular  water  (Doremus 
(46)).  The  transport  of  the  diffusing  species,  dissolved  molecular  water,  is 
again  postulated  to  involve  no  direct  reaction  with  the  network.  However,  as 
indicated  in  equation  (1),  a  strong  reaction  occurs  with  silicon-oxygen  bonds 
in  tne  network,  accounting  for  its  high  solubility  as  well  as  the  square  root 
dependence  of  the  reacted  water  concentration  on  ambient  water  vapor  pressure 
The  reaction  products  (i.e.  SiOH  groups)  are  postulated  to  be  immobile,  while 
tne  itiolecularly  dissolved  water  diffuses  with  a  diflusivity  D,  as  discussed 
above . 


09  Solubility,  Diffusivity:  Previous  Permeation  Measurements 

Previous  experimental  values  for  oxygen  uilfusivity  in  SiGp  were  deter¬ 
mined  using  variants  of  the  basic  technique  mentioned  above.  By  mass 
spectrometry,  Norton  (41)  measured  the  evolution  of  (^0)2  gas  inside 
evacuated  silica  bulbs  which  were  subjected  lor-  a  time  to  a  known  oxygen 
pressure  outside.  haul  and  Dumbgen  (51)  and  Sucov  (.52)  introduced  measured 
quantities  of  gas  into  chambers  containing  silica  fibers,  then  observed 

the  decrease  in  mole  fractions  with  a  mass  spectrometer;  Williams 

(53)  did  the  same  for  both  bulbs  and  fibers.  Their  results  (as  they  appear 
in  the  works  cited)  nave  here  been  compiled  as  an  Arrhenius  plot  in  Figure  2, 
and  are  considered  at  length  in  the  Discussion;  it  may  be  remarked  here  that 
the  three  sets  of  very  low  oxygen  Jif fusivities  were  all  measured  by  means  of 
tracer  180  (to  be  discussed  later)  while  the  other  was  measured  with  natural 
Norton  (49)  a, so  determined  the  Oi  solubility  ratio  to  be  1.9(10)~3 
(ut  ; 078  C)  ana  2.3(10)“3  (at  950  C)  and  observed  that  the  rate  of  permeation 
was  proportional  to  the  first  power  of  the  ambient  oxygen  pressure.  None  of 
Cue  other  works  presents  an  explicit  value  of  the  O7  solubility  ratio. 


46.  R,  H.  Doremus,  in  Reactivity  of  Sol  ids ,  Mitchell,  de  Vries,  Roberts 
and  Cannon,  Eds.  (Wiley,  1969)  p.  667 

49.  F.  a'.  Norton,  Nature  191  ,  701  (1961) 

50.  F..  Rosencher,  A.  Straboni,  S.  Rigo  and  G.  Amsel,  Appl  Phys  Lett  34 , 
254  (1979) 

51.  R.  Haul  and  G.  Dumbgen,  7,  Elect  rochem  66_,  636  (1962) 

52.  E.  W.  Sucov,  J  Am  Ceram  Soc  4_6,  14  (1963) 

33.  E.  L.  Williams,  .1  Am  Ceram  Soc  48,  190  (1965) 


u  r  e  2 . 


P 


(Following  Page)  Temperature  dependence  of  D 
revious  permeation  results. 


ef  f  :H20 


and  D02 


KEY  TO  FIGURE 

SYMBOL  SOURCE  Dq  (cm2/s)  Q  (kcal/mol) 


• 

Norton  (49) 

2.9(10)“4 

27.0 

— 

Haul  and  Dumbgen  (51) 
as  quoted  in  Terai  and 
Hayami  (92) 

1.2(10)~5 

56.0 

□ 

Sucov  (52) 

1.5(10)~2 

71.2 

O  A 

Williams  (53) 

2.0(10)~9 

29.0 

X 

Moulson  and  Roberts  (54) 

1.0(10)'6 

18.3 

.  .  NYirton,  Nature  191,  701  (1961) 

*  T  Dumbgen,  7.  F.l  ectrochem  66,  636  (1962) 

.  W.  Sucdv,  Am  Ceram  Soc  46,  14  (1963) 

.  Williams,  .1  Am  Ceram  Soe  48^,  IPO  (1965) 

J,  Moulson  and  J.  P.  Roberts,  Trans  Faraday  Soc  5_7,  1208  (1961) 

R.  C.  Weast,  Handbook  of  Chemistry  and  Pbysics,  55th  ed  (CRC  Press, 
1974) 
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and  DO0  in  SiO„:  previous  permeation  results 


H  ,0  Solubility,  Di  t"  fusivitv  :  Previous  Permeation  Measurements 

In  the  ease  of  water,  Moulson  and  Roberts  (54)  measured  the  effective 
k  eiflusiVLty  in  vitreous  silica  slabs  by  observing  permeation-induced 
nges  in  optical  density  of  the  OH  infrared  absorption  peak  at  2.7  um.  By 
s  technique  they  also  determined  the  number  of  hydroxyl  groups  per  Si02 
•.<>  : o  range  from  6  x  10~^  (0600  C)  to  3  x  10~^  (01200  C)  at  1  atmosphere 

•served  it  to  vary  as  the  square  root  of  the  ambient  water  vapor 

ui..  Their  solubility  resuLts  were  confirmed  by  Shackelford  et  al  (55) 

■  «  similar  technique.  Shackelford  and  Masaryk  (45)  later  analyzed  this 

i ni’  a  slightly  temperature-dependent  activation  energy  for  OH-group 
’!i  .  The  observed  pressure  dependence  of  f 0H(Sa  t)  ,  as  well  as  its  large 

•  .  ;  caused  by  the  strung  reaction  of  water  with  Si02,  as  noted. 


cp t rat  ion  of  dissolved  molecular  water  is  inaccessible  to  direct 
'  1  sc  much  sma’ler  than^^Qr(Sat) ;  ratlu-r  it  must  be 
••••do;  i;  it,  among  known  solubility  ratios  of  nonreactive  gases 
•  its  meiecdar  size.  Doremus  (17), taking  the  molecular  diameter 
is  viscosities  at  700  C)  to  be  3.3  X  (Doreums ,  P.  133  (8),  con- 
t be  solubi* ity  ratio  of  molecular  water  is  0.01.  Using  Grove's 


2. 25 (10) 22  (cm  3) 


(7) 


•to  that  the  results  of  Moulson  and  Roberts  (54),  and  Doremus  (17) 


/H20(Sat) 

(Sat) 


\j P(atm)  x 


(@  600  C) 
(0  1200  C) 


(3) 


•cremes,  Glass  Science  (Wiley,  1973) 

.  11 .  Doremus,  J  Phvs  Ch»m  80,  1773  (1976) 

:.  F.  Shackelford  and  .1.  S.  Masaryk,  .1  Noncryst  Sol  2^,  55  (1976) 

d.  Moulson  and  J.  P.  Roberts,  Trans  Faraday  Soc  57_,  1208  (1961) 

J.  F.  Shackelford,  R.  L.  Studt  and  R.  M.  Fulrath,  J  Appl  Phys  4_3,  1619 
(1972) 
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The  values  01  Deff  which  were  measured  for  H?0  diffusion  in  SiC>2  by 
Moulson  and  Roberts  (54)  appear  in  Figure  2.  If  these  are  taken  in  conjunc¬ 
tion  with  their  measured  values  for  ^QH(Sat),  values  for  the  diffusivitv  of 
molecular  water  can  be  calculated  by  use  of  equations  (5)  and  (8),  namely 

f 8 .9(10)  9  (@  600  C)  9 

D  =  /  cm  /s  .  (9) 

\5,2(10)“7  (@1200  C) 

Finally,  the  tracer  diffusivity  may  be  calculated  on  the  basis  of  free 
exchange.  Using  equation  (6)  and  the  inequalities  /si02 >>  /^O 
we  find  that  D*  is  related  to  Deff  by 

~JoH(Sat\^  f  15(10)"^  @  600  C  .  (10) 

free  exchange  4^>Si02  w.5(10)  4  @1200  C 

The  value  of  P  0 p  sjO?  ^as  keen  taken  from  the  solubility  measurements  of 
Moulson  and  Roberts. 

0.,  and  K,.0  Diffusion  through  FiO,,  Films:  Thermal  Oxidation  Kinetics 

The  diffusion  (■:  0^_  ^nu  h2°  through  Si02  is  intimately  involved  in  the 
thermal  oxidation  of  :,i .  i.\,  .  When  silicon  is  maintained  at  an  elevated 
temperature  in  an  uy.  ■  ■:  i  -J  tg  .  tmosphere  such  as  oxygen  or  steam,  a  layer  of 
S 10  .s  formed  or,  •  : ;  urfac  ,  which  ,  r<  ws  continually  in  the  presence  of  the 
o;.:.c«nt,  Varies;,  ;  rac  •  *  experiments  cited  xn  Deal  and  drove  (56)  have  estab¬ 
lished  tr.at  oxidai  io-  ,>,  ~  <„<is  by  the  inward  transport  of  oxygen  through  the 

layer  as  u  consult  -er.t  ot  molecular! y  dissolved  oxygen  or  water.  The  silicon 
substrate  is  converted  to  new  oxide  by  the  incorporation  of  atomic  oxygen  at 
the  oxide-silicon  interface  after  molecular  dissociation.  The  layer  thickens 
by  formation  of  new  oxide  underneath  the  existing  layer. 

As  discussed  previously,  the  dissociation  products  also  react  strongly 
with  the  Si02  network  only  it  the  case  of  H2O  oxidation  (46);  for  oxygen, 
the  atomic  oxygen  reacts  significantly  only  with  the  silicon  at  the  interface 
(30).  Ir.  that  case,  as  was  pointed  out  by  Blanc  (57),  an  equilibrium  exists 
at  the  interface;  atomic  oxygen  can  either  form  new  Si02  or  recombine  into 
molecular  oxygen.  Blare  was  able  to  deduce  equilibrium  constants  for  these 
reactions  so  as  to  prouuce  a  very  impressive  match  to  a  large  body  of  data  on 
thermal  growth  kinetics,  and  disprove  the  existence  of  a  hypothetical  space- 
charge  layer  which  Deal  and  Grove  (56)  introduced  to  reconcile  the  data  to 
their  linear-parabolic  model  (discussed  later). 


46.  R.  H.  Doremus,  in  Reactivity  of  Solids,  Mitchell,  de  Vries,  Roberts 
and  Cannon,  Eds.  (Wiley,  1969)  j>.  667 

50,  E.  Rosencher,  A.  Straboni,  S,  Rigo  and  G.  Amsel,  Appl  Phys  Lett  34 , 

754  (1979! 

54,  J.  Moulson  and  J.  P.  Roberts,  Trans  Faraday  Soc  5_7,  1208  (1961) 

5b.  B.  F,.  Deal  aid  A.  S.  Grove,  J  Appl  Phys  _3£,  3770  (1965) 

37.  J.  Blanc,  Appl  Phys  Lett  33,  424  (1978) 
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It  may  tie  mentioned  that  other  competing  models  ior  the  thermal  ox  i  da  Lien 
process  also  exist,  including  transport  of  charged  interstitial  oxygen  and 
holes  (58),  transport  of  0"",  0~  and  02~  ions  through  a  nonstoichiomel  r  i  <_• 
interphase  region  beneath  the  oxide  layer  (59),  and  transport  of  excess 
oxygen  and  silicon  centers  (60).  None  of  these  provides  an  acceptable  fit 
to  the  experimental  data.  A  model  of  interfacial  dissociation/recomb J nit  ion 
similar  to  Blanc's  was  proposed  by  Ghez  and  Van  der  Meulen  (61);  this 
involved  the  incorporation  of  both  atomic  and  molecular  oxygen. 

The  growth  kinetics  of  thermal  oxide  films  provide  an  independent  method 
of  measuring  oxygen  and  water  diffusion  in  SiC>2.  By  equating  the  steady- 
state  fluxes  of  oxidant  entering  the  oxide  layer,  diffusing  through  it,  and 
forming  new  oxide  at  the  interface.  Deal  and  Grove  (56)  were  able  to  sli'w 
that  the  thickness  of  the  growing  oxide  layer  obeys  a  linear-parabolic  equa¬ 
tion  of  the  form 

x  2  +  Ax  =  B(t  +  t)  (11) 

o  o 

.'here  x  is  the  thickness  of  the  layer  at  time  t.  Their  equation  (12b) 

:  ••(ting  the  parabolic  rate  constant  B  with  the  effective  diffusivity  of  r.ht- 
oxidizing  species  appears  in  our  notation  as 


BH20  =  °e 


\  “2D0 


f/70H(Sat)//?Si02 

(12) 

p  (Sat)  , p 

2/o2  y^sio2. 

(13) 

Deal  and  Grove  measured  the  parabolic  rate  constant  B  to  have  activation 
energies  of  28.3  and  16.3  keal/mol  for  O2  and  H2O  oxidation;  by  combining 
their  measurements  of  B  with  the  values  for  Deff  found  by  Norton  and 
Moulson  (49)  and  Roberts  (54)  respectively,  they  found  that^^oqj(Sat)  an(j 
p<>..  (Sat)  agreed  with  the  permeation  values  to  within  about  10%.  They  found 
<5  lo  be  linearly  proportional  to  ambient  gas  pressure  in  both  cases.  This  is 
be  expected  from  the  observed  square-root  dependence  of  both  Deff  and 

on  pressure  for  water  and  their  linear  dependence  for  oxygen,  as 
was  pointed  out  by  Dor emus  (17). 


’7,  R „  H.  Dor emus,  J  Phys  Chem  80,  1773  (1976) 
o'  v .  .!.  Norton,  Nature  191  ,  701  (.1961) 

06.  J.  Moulson  and  J.  P.  Roberts,  Trans  Faraday  Soc  57_,  1208  (1961' 

b;*.  E,  Deal  and  A.  S.  Grove,  J  Appl  Phys  36,  3770  (1965) 

3°..  E,  M.  Fowkes  and  F.  H.  Hielscher,  Electrochem  Soc  Abstract  #182, 
Spring  Meeting,  Seattle  (21  May  1978) 

59.  A.  Lora-Tomayo,  E.  Dominquez,  E.  Lora-Tamayo  and  J.  Llabres, 

Appl  Phys  17,  79  (1978) 

60.  R.  J.  Maier,  "A  Study  of  S102  Growth  Mechanismj"  AFWL-TR-76-228 
(US  Air  Force,  Kirtland  AFB,  NM,  1977) 

61.  R.  Ghez  and  Y.  J.  Van  der  Meulen,  J  Electrochem  Soc  119 ,  1100  (1972) 


Since  che  work  of  Deal  and  Grove,  a  large  number  of  other  experiments  on 
the  growth  kinetics  of  thermal  oxides  have  been  conducted.  Among  these  are 
measurements  in  oxygen  at  1  atm  (62,  65,  64 ,  66,  40,  66,  67,  68)  and  steam 
(69,  70);  the  activation  energies  reported  by  all  these  investigators  are 
in  good  agreement  with  Deal  and  Grove.  Some  of  this  data  is  displayed  in 
Figure  3,  where  it  may  be  compared  to  the  permeation  measurements  shown  in 
Figure  2. 

in  addition  to  these,  a  number  of  studies  have  been  made  of  the  depen¬ 
dence  of  B  on  tne  partial  pressure  of  ambient  water.  Ota  and  Butler  (71) 
found  a  linear  dependence  at  1230  C  for  H2O  partial  pressures  between  0.1  and 
1  atm,  using  water  vapor  either  alone  or  in  mixtures  with  H2  and  Ar.  Deal 
and  Grove  also  found  this  dependence  between  0.02  and  0.2  atm  at  1100  C  with 
O7,  and  Ehara  et  al  (68)  between  0.02  and  1  atm  at  1100  C  with  O2.  Katz  and 
Howells  (3)  found  B  to  be  twice  that  expected  at  20  atm  at  725  C  using  pure 
HoO.  Irene,  however,  (64,  40)  has  found  that  B  increases  sharply  with  the 
first  traces  of  HoO  in  O2  up  to  125  ppm  followed  by  a  smaller  rate  of 
increase  for  H-iO  concentrations  from  0.001  to  0.02  atm.  The  measurements 
performed  in  this  study  contradict  those  of  Irene  -  they  imply  a  purely 
linear  pressure  dependence  down  to  80  ppm.  A  likely  explanation  for  the 
discrepancy  and  suggestions  for  its  resolution  are  put  forth  at  length  in 
the  Discussion  section. 


3.  1..  E.  Katz  and  B.  F.  Huwelas,  J  El.-ctrochem  Soc  126,  1822  (1979) 

40.  E.  A.  Irene  and  R.  Chez,  J  Electrocbem  Soc  124 .  1757  (19771 

62.  A.  G.  Revesz  and  R.  J.  Evans,  J  Phys  Chem  Solids,  30,  551  (1969) 

63.  M.  A.  Hopper,  R.  A.  Clarke  and  L.  Young,  J  Electrochem  Soc  122,  1216 
(1975) 

64.  E.  A.  Irene,  J  Electrochem  Soc  121,  1613  (1974) 

65.  E.  A.  Irene  and  Y.  J.  Van  der  Meulen,  J  Electrochem  Soc  123,  1384  (1976) 

66.  B.  E.  Deal,  J  Electrochem  Soc  125 ,  576  (1978) 

1 ,  / .  11.  E .  Heal,  D.  W.  Hess,  .1.  I).  Plummer  and  C.  )’.  Ilo,  .1  Electrochem  Soc 

125,  339  (1978) 

68.  E.  Ehara,  K.  Sakuma  and  K.  Ohwada,  J  Electrochem  Soc  126,  2249  (1979) 

69.  7.  Nakayama  and  F.  C.  Collins,  J  Electrochem  Soc  113,  706  0966) 

70.  k.  A.  Pliskir. ,  TBM  J  Res  Dev  \0,  198  (1966) 

71.  Y.  Ota  and  S.  R.  Butler,  J  Electrochem  Soc  121 ,  1107  (1974) 
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Figure  3.  (Next  page)  Previous  measurements  of  parabolic  rate 
constants  for  thermal  oxidation  of  Si  in  O2  and  steam.  The  values 
are  those  given  by  authors  for  oxidation  at  1  atm*  Where  necessary, 
the  quoted  values  were  adjusted  to  1  atm  by  assuming  that  B  was 
linearly  proportional  to  gas  pressure. 
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j’reviou.s  r  e.n-uirement  s  of  parabol  io  rate  constants  tor  thermal 
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cs  -v.;eii.  .Note  similarity  of  activation  energy  to  nermeation  values. 


■  ,  '>:_i  1  i.jig  Techniques 

•  •  ir-aurt'aco  chemical  analysis  using  ion  beams  is  by  now  a  common  practice 
■.i-rials  research;  see,  Cor  example,  Thomas  and  Cachard  (72),  Mayer  and 
f  ( 7  3 )  and  Meyer  et  al  (74).  The  technique  most  often  employed  is 
'lord  Hackscattering  Spectrometry  (RBS)  (75),  in  which  iow-Z  ions  (at  a 
fixed  energy  somewhere  below  10  MeV)  are  directed  at  a  sample.  After 
undergo  large-angle  nuclear  (Rutherford)  scattering  at  some  depth  within 
'mpie,  tire  ions  themselves  are  detected  by  a  detector  capable  of  .inn!yz- 
•c i r  energy.  Since  the  ions  lose  energy  throughout  their  path  in  the 
'  <  ,  the  energy  which  remains  upon  their  emergence  (at  a  given  neat tor- 
■  1  e )  serves  to  indicate  the  depth  at  which  the  scattering,  occurred. 

•.•rgv  spectrum  of  backseat  ter  ed  ions  indicates  the  distribution  of 
j  c-s .  It  should  be  noted  that  since  the  elastic  scattering  cross 

is  a  siowlv  varying  function  of  energv,  the  probability  of  scattering 
e:o:ig  function  oi  depth  (for  near-surface  events).  RBS  finds  a 

ion  in  determining  distributions  of  high-Z  impurities  ir  ■<■■*■- 
i  he  'targe  difference  in  scattered  ion  energy  fol  i  owing  high-?, 
v  /  sc.  tiering  events  makes  the  two  types  ot  scattering  ro'di'v 

‘.file  the  7,  dependence  of  the  Rutherford  cross  section  supprr  s.'t’1 
>  but  !on  of  the  low-?  matrix. 

i  imin.it  ion  ol  oxygen  d  i  si  r  flint  ions  ,  however,  the  ana  !  vs  i  s  teeli- 
ed  in  this  study  was  that  of  nuclear  resonance  profiling  (/(>,  77'. 
c  ique,  as  applied  to  !  ^0  profiling,  is  illustrated  in  Figure  4,  A 
’ '-meed  beam  o!  protons  from  a  Van  de  Graaf  electrostatic  accelerator 
oxiue  layer  after  being  energy  analyzed.  The  cross  section  for 
■•"'if  reaction  ^n(p ,  of  )  1  ’N  between  the  protons  and  the  ^0,  namely 

18  IS 

P+  0  -  r+o<4-  3.980  MeV 

iPaunced  narrow  resonance  et  a  proton  energy  of  829  keV  (Mayer  and 
(73))  .  Protons  enter  at  a  higher  energv  than  that,  and  lose 


' ,  (homos  and  A.  Cachard,  Material  Characterization  I'sing  Ion  Beams 
■■vs,  19  78) 
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Figure  4.  Principle  of  nuclear  resonance  profiling.  The  graph  shows  the 
cross  section  vs  incident  proton  energy  at  8lab  =  150°  for  180(p,O<)15N 
reaction.  The  peak  energy  is  at  629  keV;  protons  enter  with  slightly  higher 
energy.  The  sensitive  region  is  located  at  the  depth  where  the  proton  energy 
has  fallen  to  that  value.  Alpha  particles,  emitted  with  an  energy  of 
3.4  MeV,  easily  escape.  (From  Mayer  and  Rimini,  p  163  (73)) 


energy  continuously  along  their  path  in  the  layer  until  they  reach  the 
resonance  energy,  at  which  depth  the  reaction  occurs.  (The  reaction  can,  of 
course,  occur  at  any  other  depth,  but  with  much  less  likelihood.)  Alpha 
particles  produced  via  the  reaction  escape  from  the  surface  and  are  detected; 
their  yield  (relative  to  the  proton  fluence)  indicates  the  concentration 
at  that  depth.  By  systematically  varying  the  energy  of  the  proton  beam  and 
observing  the  alpha  particle  yield,  the  concentration  profile  can  be 
determined. 
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In  order  to  determine  an  unknown  0  concentration  profile  by  the  tech¬ 
nique  of  nuclear  resonance  profiling,  one  observes  the  yield  of  alpha  parti¬ 
cles  (per  steradian  per  incident  proton)  from  the  reaction  of  protons  with 
’-»0  distributed  within  the  oxide  layer.  A  general  expression  for  the  alpha 
particle  yield  (vs.  proton  beam  energy)  is 


(Eb)  = J dxjfl(x)  J  & E£g(Eb,  Ei) J dE^(E)  f(E,  Ei,  x)  (14) 


where 

E^  is  the  mean  energy  of  the  incident  proton  beam, 

E  is  the  energy  of  the  reacting  protons, 

E.j  is  the  actual  energy  of  protons  incident  on  the  surface  of  the  layer, 

y  is  the  thickness  of  the  oxide  layer, 

x  !s  the  depth  at  which  the  reaction  occurs, 

18 

is  the  0  concentration  profile  (i.e.  atoms  per  unit  volume  at  depth  x), 

c—  is  the  differential  cross  section  of  the  1^0(p,O(  )15N  reaction, 

d‘>' 

•.(Ei,,  Ei)dEi  is  the  fraction  of  protons  in  the  beam  which  have  an  incident 
•lergy  between  Ef  and  E±  +  dE^, 

(E,  Ei,  x)dE  is  the  probability  that  a  proton  which  has  started  with  energy 
j  and  arrived  at  depth  x  will  have  an  energy  between  E  and  E  +  dE. 

If  the  functions  g(Eb,  E^) ,  f(E,  E^,  x)  and  -g^-(E)  are  all  previously  known 
and  the  yield  Y(Eb)  is  observed,  then  the  function^(x)  can  in  principle  be 
determined  by  inversion  of  the  integral  equation. 
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An  enlightening  discussion  of  this  matter  appears  in  Dunning  et  al  (78),  in 
which  a  similar  calculation  (involving  the  27Al(p,  y  )2®Si  reaction)  appears. 

It  should  be  noted  that  the  energy  of  the  emerging  alpha  particles  is 
(a)  considerably  greater  than  that  of  the  incident  protons,  and  (b)  almost 
equal  to  the  Q-value  of  this  inelastic  reaction.  They  lose  at  most  about 
40  keV  (having  been  emitted  with  2.4  MeV)  in  traversing  the  2000  &  thick 
Si02  layers  in  this  study  (Mayer  and  Rimini,  p.  16  (73)). 

Applications  of  nuclear  resonance  profiling  of  ^®0  have  included  oxygen 
diffusion  in  Zr02  (79),  silicon  (80),  GaP  (81),  Ti02  (82;,  HO2  and  Cr203 
(83),  and  water  diffusion  in  Ta  oxide  (84). 

Nuclear  resonance  profiling  contrasts  with  RBS  in  several  important 
respects.  First,  since  the  depth  at  which  the  reaction  occurs  is  determined 
solely  by  the  incident  particles'  energy,  the  detection  system  need  not  pro¬ 
vide  energy  analysis  for  the  reaction  product  (except  to  screen  out  products 
of  competing  reactions).  Second,  the  specificity  of  the  reaction  allows  the 
detection  of  minute  quantities  of  an  isotope  in  the  presence  of  any  other 
material  (subject  to  the  non-interference  of  other  reactions).  This  makes  it 
especially  valuable  for  detection  of  low-Z  materials,  whose  resonant  reac¬ 
tions  tend  to  be  readily  distinguishable.  Third,  and  most  important,  narrow 
resonances  provide  high  resolution  in  depth  -  the  2.1  keV  width  of  the  629 
keV  I80  resonance  corresponds  to  a  depth  spread  of  350  &  (as  is  discussed  in 
greater  detail  later) . 
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It  is  important  to  distinguish  nuclear  resonance  profiling  from  the  other 
principal  profiling  technique  based  on  direct  observation  of  nuclear  reac¬ 
tions  (85,  86).  That  technique,  pioneered  by  Amsel  (1963),  uses  reactions 
in  an  energy  domain  where  the  reaction  cross  section  is  a  slowly  varying 
function  of  energy.  As  in  RBE,  the  energies  of  the  emerging  reaction 
products  are  analyzed;  the  depth  distribution  of  reacting  material  is 
provided  by  the  resulting  energy  spectrum.  This  technique  has  been  employed 
in  a  number  of  studies  involving  oxygen,  including  1^0  diffusion  in  Alpha- 
zirconium  (87,  88),  170  diffusion  in  zirconium  oxide  (89),  1®0  diffusion  in 
quartz  (90),  and  determination  of  1^0  concentration  profiles  in  Si02  films 


(91). 
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EXPERIMENTAL  PROCEDURE 


Materials  Preparation 

Single-crystal  silicon  wafers  (2"  diam)  were  obtained  from  Monsanto 
Corp.  They  were  of  (100)  orientation,  with  one  face  polished;  their  resis¬ 
tivity  was  4  ohm-cm  (n-type) .  They  were  first  degreased  in  boiling  trichlor- 
ethylene,  cleaned  in  both  H2SO4  (at  90  C)  and  HF,  and  rinsed  in  deionized 
water.  They  were  then  oxidized  to  2000  A  in  flowing  dry  oxygen  at  1  atm  at 
1200  C  in  a  commercial  3-zone  tubular  furnace  (Thermco)  which  was  lined  with 
a  quartz  diffusion  tube.  Both  the  sample  preparation  and  oxidation  proce¬ 
dures  are  standard  in  the  semiconductor  processing  industry. 

1 8 

Following  the  oxidation,  0  was  implanted  into  the  Si02  films.  The 
*®0  was  in  the  form  of  isotopically  enriched  gas  (70  atom  %,  obtained  from 
ProChem  Isotopes) ,  a  lecture  cylinder  of  which  was  mounted  in  a  commercial 
ion  accelerator  (at  RCA  Laboratories,  Princeton,  N.  J.).  This  device  ionized 
the  gas,  then  electrostatically  accelerated  the  ions  to  the  desired  energy 
through  a  mass  analyzer  which  accepted  only  ions  of  the  desired  atomic  weight 
The  implantation  was  accomplished  by  raster-scanning  the  ion  beam  across  the 
polished  faces  of  the  wafers  to  provide  a  uniform  fluence  of  3(10) ^ ^  ions/cm^ 
To  avoid  H2^0  contamination,  the  mass  analyzer  was  set  at  36  AMU  rather  than 
18  AMU,  so  that  the  projectile  ions  were  actually  (^®0-  0)+.  To  reach  the 

desired  depth,  the  molecular  energy  was  set  at  80  keV  per  nucleus).  A  mass 
analysis  of  the  residual  gas  within  the  ion  implantation  apparatus^  which  ap¬ 
pears  as  Figure  5,  showed  the  major  contaminants  to  be  (^0-  °0),  ^2  and 
^uAr,  in  quantities  insufficient  to  affect  the  fluence  measurement  signifi¬ 
cantly.  The  wafers  received  no  immediate  post-implant  treatment. 

Finally,  samples  were  prepared  by  cleaving  the  oxidized  and  implanted 
wafers  into  5mm  x  10mm  sections  which  had  been  individually  labeled  to  indi¬ 
cate  their  original  locations  on  the  wafers. 

Thermal  Treatments 


The  general  procedure  involved  heating  samples  in  atmospheres  containing 
various  water  vapor  concentrations.  In  each  atmosphere,  samples  were  heated 
at  a  number  of  temperatures  in  order  to  determine  the  activation  energy  of 
diffusion  in  equation  (3).  Comparison  of  the  pre-exponential  factors  Do*  ob¬ 
tained  in  different  atmospheres  at  a  single  temperature  provided  the  means  of 
determining  the  pressure  dependence  of  D* .  The  extent  of  diffusion  as  a 

function  of  heating  time  at  a  fixed  temperature  and  atmosphere  was  observed 
in  order  to  determine  reaction  rates  involved  in  water  permeation. 

Temperature  Dependence 

The  first  set  of  thermal  treatments  was  carried  out  in  room  air.  .The 
room  was  maintained  at  a  temperature  of  21  C  and  a  relative  humidity  of  ap¬ 
proximately  50%.  Samples  were  placed  in  a  horizontal  3-zone  tubular  alloying 
furnace,  (Pacesetter  II,  Thermco  Products  Corp.)  lined  with  a  6.4  cm  ID  X  91 
cm  long  fused  quartz  diffusion  tube,  open  at  one  end.  The  furnace  was 
equipped  with  a  proportional  temperature  controller  (Ana-Lock  Series  321, 
Thermco  Products  Corp.)  for  each  of  its  three  integral  Type  K  (Ni-Cr /Ni-Al) 
thermocouples.  Sample  temperatures  were  monitored  by  means  of  a  separate 
type  R  (Pt/Pt-1 3%Rh)  thermocouple  which  was  inserted  in  the  diffusion  tube. 
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The  thermocouple  wires  were  fixed  within  a  91  cm  long  X  0.6  cm  OD  quartz 
tube;  the  hot  junction  was  located  near  its  closed  end.  The  thermocouple 
indicator  (Model  DS-350,  Doric  Scientific  Corp.)  had  been  factory-calibrated 
for  a  type  R  thermocouple;  it  had  a  digital  display  with  a  resolution  of 
0.1  C.  The  furnace  controls  were  adjusted  while  monitoring  the  temperature 
at  various  points  so  that  an  8  cm  long  "flat-zone"  was  obtained  near  the 
center.  Diffusions  were  carried  out  by  loading  a  batch  of  samples  on  a  4.1 
cm  X  15  cm  plate-type  slotted  carrier  made  of  fused  quartz  (obtained  from 
Quartz  International  Corp.)  which  had  been  preheated  for  at  least  45  minutes 
at  the  center  of  the  furnace.  The  carrier  was  pulled  out  to  the  mouth  of  the 
diffusion  tube  using  a  hooked  quartz  rod,  the  samples  were  loaded  quickly  and 
the  carrier  returned  to  a  position  such  that  the  samples  were  within  the  flat 
zone;  the  thermocouple  was  then  inserted  adjacent  to  the  samples.  Sample 
placement  and/or  removal  involved  withdrawing  the  carrier  from  the  furnace 
for  no  more  than  half  a  minute  per  operation.  Samples  were  transferred  to  a 
room-temperature  pyrex  surface  where  they  cooled  rapidly.  Samples  were  al¬ 
ways  placed  with  the  implanted  face  up.  Timing  was  done  with  a  quartz-cry¬ 
stal  digital  wrist  watch. 

The  next  set  of  thermal  treatments  was  carried  out  in  dry  nitrogen  at 
1  atm.  The  same  furnace  controller,  thermocouple  and  diffusion  tube  were 
used.  A  fused  quartz  diffusion-type  end  cap  was  placed  on  the  gas  supply 
end  of  the  diffusion  tube;  the  other  end  was  left  open.  Gas  was  obtained  via 
boiloff  of  liquid  nitrogen;  flow  was  maintained  using  a  tube  flowmeter  (Sho- 
Rate,  Brooks  Instrument  Company)  at  500  cc/min  (equivalent  at  70°  F,  1  atm). 

A  similar  procedure  was  followed  for  thermal  treatments  in  steam  at 
1  atm.  in  this  case,  however,  the  furnace  at  first  was  a  horizontal  single¬ 
zone  tubular  resistance  furnace  (Model  54233,  Lindberg/General  Signal),  lined 
with  a  5.1  cm  ID  x  91  cm  long  fused  quartz  diffusion  tube.  The  furnace  was 
equipped  with  an  integral  proportional  temperature  controller  (Model  59544, 
Lindberg/General  Signal)  for  the  furnace's  integral  thermocouple.  After  a 
time,  the  furnace  and  controller  were  replaced  with  a  horizontal  single-zone 
tubular  resistance  furnace  (Model  1127,  Marshall  Furnace  Co.)  lined  with  a 
3.8  cm  ID  x  117  cm  long  fused  quartz  diffusion  tube.  Furnace  temperature  was 
controlled  using  a  proportional  controller  (Eurotherm  Model  901-2075,  Marshall 
Furnace  Co.)  and  a  type  K  taermocouple  whose  junction  was  located  at  the  cen¬ 
ter  of  a  small  tube  bored  parallel  to  the  furnace  tube.  Sample  temperatures 
were  monitored  by  moans  of  a  separate  type  K  thermocouple  whose  hot  junction 
w.is  located  near  the  closeU  end  of  a  91  cm  long  x  0.95  cm  OD  quartz  tube, 
which  was  inserted  in  the  diffusion  tube.  The  thermocouple  EMF  was  read  out 
on  the  200  mV  scale  of  a  3-1/2  digit  multimeter  (Model  8020A,  John  Fluke 
Corp.),  and  converted  to  temperature  with  the  aid  of  appropriate  tables 
(Omega,  1979).  The  readings  were  compensated  for  cold  junetio.  temperature 
by  means  of  an  electronic  cold-junction  compensator  (Model  MCJ-K,  Omega  En¬ 
gineering  Corp,),  The  type  K  and  R  thermocouples  were  c ross-checked  by  put¬ 
ting  both  into  the  furnace  at  various  temperatures  between  550  C  and  600  C. 

The  type  K  thermocouple  was  calibrated  by  separating  the  two  wire  ends 
and  bridging  them  with  a  small  horizontal  strip  of  silver  foil  whose  surfaces 
had  been  cleaned  by  abrasion.  The  ends  of  the  silver  strip  were  crimped  onto 
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the  wires,  which  were  then  heated  under  high  vacuum.  The  last  DMM  reading  be¬ 
fore  the  strip  melted  corresponded  to  a  temperature  of  958  C,  which  is  with¬ 
in  4  C  of  the  melting  point  of  silver  (1).  This  discrepancy  was  judged  to  be 
within  the  error  of  reading  on  the  digital  multimeter,  especially  since  the 

temperature  was  still  increasing  at  the  time  of  the  last  reading,  so  that  the 
meter  readings  were  taken  as  indicating  the  true  temperature. 

The  steam  was  generated  by  boiling  deionized  water  in  a  stoppered  2-liter 
foiling  flask,  as  shown  in  Figure  6.  The  boiling  rate  was  governed  by  a  vari¬ 
able  autotransformer  supplying  power  to  a  hemispherical  heating  mantle  around 
t  •  flask.  Steam  was  conducted  to  the  diffusion  tube  through  a  horizontal  py- 
■x  Lube  emerging  from  the  neck  of  the  flask;  this  was  joined  to  the  narrowed 
e  oi  the  diffusion  tube  via  a  tapered  ground  joint.  The  entire  length  of 
h"  ■:  onta.l  tubing  and  the  neck  of  the  flask  were  wrapped  in  fiberglass  heat- 

and  maintained  at  a  temperature  of  120  C  to  avoid  condensation  there, 
'perature  was  monitored  with  a  Type  T  (copper-constantan)  thermo- 
•  '  .  An  .anti-hackstreaming  plug,  formed  by  a  wad  of  fiberglass  tape,  was 

the  thermocouple  tube.  When  inserted  into  the  diffusion  tube,  it 
:  ..  slight  positive  pressure  to  be  maintained  in  the  furnace.  Sample 
i  :  i.'ii,  a  procedure  which  took  approximately  15  seconds,  was  delayed  until 
.  :  >'  had  been  replaced  by  steam  and  the  furnace  temperature  had  stabilized. 

■I'les  rested  on  a  2  cm  x  6.5cm  quartz  boat  which  was  inserted  to  the 
•  of  the  furnace  with  a  hooked  steel  rod. 

Pressure  Dependence 

'i.aai.Ly,  a  set  of  thermal  treatments  was  carried  out  in  pure  water  vapor 
•  iu.cd  pressures.  These  were  performed  within  the  46  cm  diameter  x  79  cm 
hi;,;  >el i  jar  of  a  vacuum  evaporation  station  (Model  VS-400,  Veeco  Corp,),  as 
>  ■■■•;  in  Figure  7.  The  samples  were  placed  inside  a  17.6  cm  long  section  of 

.  'Jb  quartz  tubing.  A  double-hole  ceramic  rod  with  a  type  K  thermocouple 

■*..  at  one  end  was  also  inserted  in  the  tubing  so  that  the  junction  was 
•ample  location.  The  tubing  was  then  inserted  into  a  1.35  cm  ID  x  10.5 
1  tubular  resistance  furnace  so  as  to  put  the  sample  at  its  center.  The 
mount,  a  small  pedestal,  rested  on  the  61  cm  diam  stainless  steel 
o  of  the  evaporation  station.  Electric  power  for  the  furnace  was  sup- 
>  variable,  autotransformer  through  feedthroughs  in  the  baseplate; 

.  SOW  was  required  to  maintain  a  sample  temperature  of  800  C.  The  ther- 
■’e  wires  were  likewise  connected  to  electrical  feedthrough  in  the  base- 
p  .  Water  vapor  was  supplied  by  deionized  water  in  a  stoppered  vacuum 

i vsting  on  a  thermostatic  hotplate  to  maintain  a  constant  water  temper- 
The  vapor  was  conducted  through  a  length  of  copper  tubing  to  a  varia- 
■<’.  valve  (micrometer  type),  and  entered  the  vacuum  chamber  through  the 
'  ; .unce  valve.  A  continuous  flow  was  maintained  by  exhausting  the 
' hrough  a  choked-down  roughing  valve  to  a  mechanical  forepump  (Model 
1  ‘•’i,  <  •-gent -Welch  Scientific  Co.).  Chamber  pressure  was  monitored  by  a  ca¬ 
l''  ■  >  •-■  manometer  (Baratron  Pressure  Meter  Type  220-2A3-10,  MKS  Instruments, 

Ip.-  whose  sensing  tube  entered  the  chamber  through  a  port  in  the  baseplate. 
Si., via  is  from  both  the  manometer  and  the  thermocouple  were  read  via  digital 
multimeters  and  recorded  on  strip  chart  recorders  (Model  680,  Hewlett-Packard 
Corp . ) . 


1.  R.  C.  Weast,  Handbook  of  Chemistry  and  Physics,  55th  ed  (CRC  Press,  1974) 
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Figure  6.  Apparatus  for  thermal  treatments  in  steam 


thermal  treatments  at  low  pressure. 


The  thermal  treatments  were  carried  out  hv  hirst  loading  the  samples  and 
thermocouple  in  the  furnace.  The  bell  ,ar  was  then  lowered  and  the  system 
evacuated  by  the  forepump  through  the  diffusion  pump  backing  line.  For  this 
operation,  the.  high-vacuum  valve,  the  backing  valve  and  the  air  admittance 
valve  were  ail  fully  open,  while  the  roughing  valve,  the  mechanical  pump  vent¬ 
ing  valve  and  the  variable  leak  valve  were  all  shut.  The  latter  was  then 
opened  to  evacuate  the  air  from  the  water  supply  flask  and  tube  and  to  degas 
the  water,  inert  shut  to  re-evacuate  the  system.  Next,  water  vapor  was  admit¬ 
ted  :.'v  adjusting  tne  variable  leak  valve  until  the  observed  rate  of  pressure 
'  n c r>.  isc  was  such  that  the  pressure  reached  its  desired  value  in  5  min.  At 
tr  point  : no  roughing  valve  was  "cracked"  open  to  maintain  that  value  stablv. 
.*  v  is  again  evacuated,  this  time  through  the  hacking  valve,  and  the 
‘■;t  e  was  wrought  up  to  temperature.  No  difficulty  was  encountered  in  keep- 
1  i:‘  temperature  constant  to  within  about  +  3C  by  manually  ad  lusting  the 

,r  >•(  rails;  ermer  voltage.  The  pressure  was  then  brought  up  to  its  present 
‘'•'..i  througii  shutoff  o’  the  backing  valve.  On  short  runs,  pressure  excur- 
.sbius  Were  prevented  through  manual  adjustment  of  the  variable  leak  valve;  on 
■  •  can:,  thev  amounted  to  something  like  +20F  of  the  average  value. 

\U  •  1  e  ..  Re  .a;  :  I  c_  1 '  1  i  >  1  i  i  j  If  g 

c'c"  t  1  nnt.err.ine  t he  extent  of  *  0  diffusion  induced  bv  the  thermal 
'  r*  •  i  ■,  •  ’  ,u>n  particle  yield  curves  were  taken  using  the  method  of  nuclear 

i  *  a.,  vi  .sl  v  described,  in  the  Introduction.  A  proton  beam 

was  •,  ’  ae\  Van  deOraaf f  electrostatic  accelerator  (Model  A,  High 

Voltage  ;  er  mil  dorp)  .  A  portion  of  the  accelerator,  as  well  as  the  an- 
•  iii..  '  -inP.ir  '  ->.  is  shown  in  Figures  8  and  9;  prominent  features  are  indi¬ 
cate,,  in  :  ;.e  ;,r  •ceding  kev. 

•  rm  r  a  c.v  was  controlled  bv  adjusting  the  field  strength  of  a 
“I  '  .  c'i  v.  -  •  i  c.  ,<  •  magnet.  The  current  to  the  magnet  coils  was  supplied 

e,  regulated  nower  supply.  The  magnetic  field  strength  was 
p..-;cs  of  an  NMK  mussneter  (Model  M2,  Nuclear  Magnetics  Corp .  / 

,,v-  .  i  whose  probe  element  was  fixed  between  the  flat  faces  of 

■  '  Vie  gaussmete r  modulated  tne  local  magnetic  field  at  the 

•in  a  pair  oi  small  Helmholtz  coils.  The  RF  oscillator  was  set 
•  ■  .  corresponding  to  the  desired  beam  energy,  and  the  beam  energy 

.  . -a.  .  i  ;  erv’.n,.  the  magnet  current  until  the  NMR  signal  dip  centered 

/••;  *  •••  t:.e  modulation  current.  An  oscillator  frequency  of  12.050  MHz 
•n  a  Hewlf  1 1 -Packard  Model  5243L  digital  frequency  counter-)  cor- 
o  •  a  beam  energy  of  830  keV.  Excursions  in  the  magnetic  field  were 

mu  uma i .  e nougn  to  limit  beam  energy  variations  to  under  +40  eV  dur- 
.  c  ■  '■  i  i  .  min  duration  alpha  particle  count. 
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KEY  TO  PHOTOGRAPHS  OF  NUCLEAR  PROFILING  APPARATUS  (pp  33,  34) 


VAN  DE  GRAAFE  VERTICAL  DRIFT  TUBE:  Protons  travel  downward. 

FRAME  OP  ANALYZ INC  MAGNET  (Coils  are  in  dark  area  below):  Protons 
undergo  90°  deflection  with  38  cm  radius  of  curvature,  then  emerge 
■"ward  right. 

SCATTERING  CHAMBER:  D-shaped,  multiport.  On  its  rear  plate  is  a 
■rizontallv  translatable  stage  upon  whose  rotatable  center  is 

■  "ted  ‘he  sample  holder.  Protons  enter  a  front  port,  housing 

».'  ’  I  me  tor,  which  is  attached  to  the  horizontal  beam  tube. 

,  >■  ■  ->rts  *’oid  vacuum  gauge,  electrical  feedthroughs,  vacuum 
.  n  ■  system  and  viewing  window.)  Extending  through  the  top 
*'.■>.  •  :  •:  a  rotatable  f eedthrough  for  changing  detector  angle 
t  ce  outside.  Bevond  it  stands  a  swi  tchable  microammeter, 
tu  the  collimator  sectors  for  beam  alignment . 

•Vtr.E  SENSITIVE  PREAMP:  Mounted  on  rotatable  feedthrough  holding 
•  cc  tor . 

' NT  ILLATION  DETECTOR.  Nal(Tl)  crystal /photomul tipi ier  tube 
■-cmblv  and  plug-in  preamp,  shielded  with  lead  bricks,  for  gamma 
■v  inlibratton  of  beam  energy.  It  should  be  noted  that  the 

■  t  ter  "F,"  in  the  upper  right  corner  is  a  wall  plaque  denoting 
>.st"  rather  than  a  label  for  this  diagram. 

•; ! M  BIN:  Spectroscopy  and  DDL  amplifiers.  Below  this  is  the 
■  V  newer  supply  for  the  photomultiplier. 

■  i '  i  <.r  ;V''OR  UNIT:  For  NMR  gaussmeter,  connected  by  a  short  cable 

nrobe  mounted  in  magnet  gap  on  far  side.  Immediately  to  its 
■  is  a  digital  counter  to  facilitate  setting  frequency;  to  its 
ignt  is  an  RE  amplifier  for  oscillator  signal. 

ILLOSCOPE:  Used  to  display  NMR  signal.  To  its  left  (not  shown) 
magnet  current  controller,  NMR  power  supply  (modified  to  include 
■■ur, actor  for  remote  trimming  of  oscillator  frequency)  and  digital 
on-nev  counter. 

"  '  \"V'R'S  CONSOLE  FOR  VAN  DE  GRAAFF  ACCELERATOR. 

’  'CHANNEL  "ELSE  HEIGHT  ANALV7.ER . 

TARGET  CURRENT  INTEGRATOR. 
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prater:  experiment. 


Immediately  beyond  the  analyzing  magnet  was  a  pair  of  control  slits. 
These  served  to  allow  passage  of  only  those  protons  with  the  set  energy,  as 
well  as  to  provide  a  feedback  signal  for  closed-loop  fine  control  of  beam 
energy.  Beyond  these  lay  the  horizontal  beam  tube,  at  the  end  of  which  was 
a  stainless  steel  multiport  scattering  chamber,  shown  schematics! 1 v  in  Figure 
]0.  This  housed  a  sample  holder  and  a  particle  detector.  The  sample  holder 
was  a  4 .  9  cm  diameter  stainless  steel  disc  with  a  shallow  groove  milled  ■>- 
cross  its  face;  batches  of  samples  stood  side  bv  side  in  the  groove  and  were 
In- :  d  in  place  bv  small  clips  extending  over  its  edges.  F.acb  batch  included  a 
-  er'iire  sample  that  had  been  ^0- implanted  but  not  thermal  lv  treated.  The 
,  bolder  was  mounted  with  its  face  normal  to  the  proton  beam  and  biased 


V  to  co 


1 1  j 


ect  secondary  electrons;  it  was  movable 


the  detector  mount.  The  detector  was  a  silicon  surface 
device  {■'-">  m~>-  area,  100  Ura  depletion  depth,  Ortec)  located  17  mm 

at  an  angle  of  J 10°  to  the  beam  direction.  Tt  was  shielded  bv 


a  s  wa  s 
> ! 


samp  1 e 


-2 


from  outside  the 

bar¬ 
ren 
a 


.st.er  film  whose  thickness  of  1.1  mg-cm  ^  was  sufficient  to  prevent  anv 
ac k.soattereb  protons  from  reaching  it  while  allowing  penetration  of  the 


•V  aipiia 
:‘<0  (.bar 


\  •'<:  . 
reached 


lei 


lit  i  ’ 


U"i 


tiaiiJ 


The  detector  pulses,  after  being  amplified  (by  a  Moot 
Preamp  and  Model  1417  Spectroscopy  Amplifier,  both 
.  ,  w.-rv-  displayed  on  a  2048-channel  pulse  height  analyzer 
racnr-Nor t hern  Corp. ) ,  which  performed  a  "region-of-interest" 
it  -iisp  lave.;  a  running  total  of  pul  ses  accumulated  in  the  alpha- 
fit  i  :i  .»•>.  cozen  by  a  signal  from  a  current  integrator  (Model 
incicating  that  the  proton  fluence  on  the  sample  had 
set  value  (tvpicallv  50  uC).  The  "region  of  interest"  was  an 
c  . (idling  the  alpha  particle  energy;  it  was  set  on  the  PHA  so  as 
iidi.  all  the  alpha  particles  registered  in  the  narrow,  well-defined 
pe.,k.  Background  counts  were  negligible.  Alpha  particle  counts 
rc.iT.  1060  to  4000  per  50  uC  run;  beam  current  was  held  to  0.5  uA  to 
e i : i Mg  the  samples  (see  Discussion  for  clarification  of  this  point). 

t-  -am  enorgv  setting,  the  sample  holder  was  translated  to  place  each 
sun  pi  os  consecutively  in  the  1  mm-diameter  beam.  The  range  of  beam 
settings  spanned  was  typically  13  keV  on  either  side  of  the  energv  at 
1,.-  I  vie  Id  was  maximum.  This  corresponded  to  a  projected  range 
: .  •  ,• ; t her  side  at  the  location  of  the  distribution  peak  -  enough  to 
>:, tng  of  fate  entire  2000  A  thickness  of  the  oxide  layer. 


;  :  a  ;  j  g  v  a]  p'.a  vie1  d  curves  were  taken,  a  calibration  was  performed  on 

..  ae  Oraat  f  acce  lerator  in  order  to  functionally  relate  the  proton  beam 
ip  the  observed  MMR  oscillator  frequency.  This  was  done  also  by  means 
i.  i  ear  resan.n.ce  profiling,  but  with  (p ,  V>  rather  than  (p,o()  reactions, 
e  ’  i i s  the  reactions  used.)  A  thick  target  containing  a  uniform  con- 
.  .  .  ! rat  ion  of  appropriate  nuclide  was  placed  in  the  proton  beam,  whose  energv 
,;is  systematically  varied.  The  yield  of  gamma  rays  was  plotted  as  a  function 
. :.Mk  frequency.  As  the  NMF  frequency  (and  thus  the  proton  energv)  was  svs- 
t  i  ca  I  I  v  increased,  the  yield  increased  relatively  si  owl  v  as  long  as  the 
o rot  on  energv  remained  below  the  relevant  resonant  energv.  (These  gamma 
ruvs  c.-i;.ic  from  a  mu  i  t  i  tude  of  smali  resonances  in  the  (p,)()  cross  section.) 
'••wevor ,  is  -non  as  the  proton  energy  reached  the  main  resonant  energy,  the 
■  ■  i  pc  r  i  .a  sed  sudden  I  v ,  reflecting  the  increased  number  of  interactions  at 
1  t  •  .  ,  .  v.  As  the  proton  energy  was  further  increased,  the  yield  resumed 
.  isc.  retirci  inn  the  slow  variation  in  the  reaction  cross  section, 
example  of  a  t  hick-target  gamma-ray  vield  curve  is  shown  in  Figure  11. 
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;ure  10.  Nuclear  scattering  chamber.  Successive  samples  are  irradiated 
horizontally  translating  the  back.  Target  current  is  collected  by  a 
ip  ring  to  allow  rotation.  See  KEY  TO  PHOTOGRAPHS  for  more  description. 


TABLE  1-  ACCELERATOR  ENERGY  CALIBRATION  REACTIONS 


Observed 


Rea 

ct  ion 

Ep  (keV) 

T  (keV) 

E^  (MeV) 

C T  (mb) 

V 

— nmr 

(MH: 

19,,, 

r  (p, 

octf)16  0 

340.5 

2.4 

6.13 

160 

8. 

930 

7  .  , 

Li  (p , 

V)8Be 

441.4 

12.2 

17.64 

6 

10. 

230 

19  . 
i  (p, 

<T> 

o 

872.1 

4.7 

6.13,  6.92 

540 

14. 

200 

27  A1  , 
Ai  (p 

,  V )28Si 

991.9 

0.10 

12.54 

- 

15. 

086 

Reference:  Mayer  and  Rimini  (2) 


f.  W,  Mayer  and  E.  Rimini,  Eds.,  Ion  Beam  Handbook  for  Materials  Analysis 
(Academic,  1977),  pp  207,  281 

37 


The  calibration  proceeded  as  follows:  A  target  thick  compared  to  a 
1  MeV  proton  range  was  prepared.  In  the  case  of  A1  this  was  a  small  block  of 
1/8"  stock;  in  the  cases  of  Li  and  F  it  was  a  layer  of  LiF  formed  by  fusing 
LiF  crystals  on  a  copper  disk.  The  target  was  placed  on  the  sample  holder 
arid  a  scintillation  detector  (a  5"x5"  Nal(Tl)  crystal  attached  to  a  5"  photo¬ 
multiplier  tube)  was  placed  on  the  scattering  chamber  nearby.  Lead  bricks 
were  positioned  upstream  to  shield  the  crystal  from  gamma  rays  coming  from 
various  beam  stops.  As  before,  an  NMR  frequency  was  selected,  and  detector 
pulses  were  displayed  on  a  multichannel  pulse  height  analyzer,  where  those  in 
an  appropriate  energy  range  were  accumulated  until  a  preset  target  charge  had 
been  reached.  The  appropriate  gamma  ray  energy  range  was  selected  on  the 
basis  of  a  detector-analyzer  system  calibration  (using  ^Co  and  *-^Cs  stan- 
d ird  reference  sources)  to  include  the  principal  gamma  rays  emitted  in  the 
reactions  listed  in  Table  I.  The  resulting  thick-target  gamma  ray  yield  was 
,  and  the  accelerator  energy  was  calibrated  by  Identification  of  the 
■  n.ency  at  which  the  step  in  the  yield  curve  occurred  (i.e.,  the  last 
-  ■  '  '  ri  Table  T)  as  the  proton  energy  corresponding  to  the  resonance. 

■■  ince  XYR  frequency  is  proportional  to  the  analyzing  magnet's  field 
cgth  and  the  magnetic  force  is  equal  to  the  centrifugal  force  on  the 
•  ••tons*  the  proton  energy  is  proportional  to  the  square  of  the  NMR  frequency. 
•  'nunc!  to  be  approximately 

E,  (keV)  =  4.36  v2  (MH  ) . 
b  z 

'Tin  cause  of  inaccuracy  is  the  nonuniformity  of  the  field  over  the  di- 
.  is  of  the  iron  pole  pieces:  the  probe,  being  of  necessity  located  near 
.  •.  1  i  inge ,  could  not  accurately  reflect  conditions  along  the  beam  path.  It 
ST"’1 1 d  be  stressed,  however,  that  this  experiment  is  self-calibrating  -  the 

■ . .  width  of  the  1H0(p,of ) reaction  allows  the  reference  sample  (im- 

r  1  -luted  but  unheated)  to  provide  an  energy  calibration  for  every  other  sample 
•  •••■  chamber. 

The  beam  was  not  purely  monoenergic:  due  to  differences  in  trajectories 
was  a  distribution  of  incident  proton  energies.  This  was  determined 
:•  "  hi-  gamma  ray  yield  calibration  curve  for  the  ~  A1  (p,y)  ^Si  reaction 
'ii  ’•  igtire  11.  Noticeably,  the  step  is  not  sharp,  but  has  a 
1l‘-  width  (here  on  the  order  of  1  keV) .  Since  the  reaction  width  is  less 
?  ha"  0.1  keV,  this  deviation  is  almost  entirely  attributable  to  the  energy 
b  of  the  beam. 
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The  thermal  treatments  given  to  all  the  processed  samples,  and  the  nu¬ 
clear  resonance  profiling  data  reflecting  the  resultant  changes  in  the 
distributions,  are  summarized  in  Table  II  and  Figures  12a  -  12c.  Table  II 
lists  the  thermal  treatments;  for  each  sample  it  shows  the  temperature,  dura¬ 
tion  and  type  of  atmosphere  employed.  It  also  shows  some  nuclear  profiling 
results  -  the  yield  of  alpha  particles  (per  50  uC  of  incident  protons)  re- 
corded  at  the  proton  beam  energy  for  which  this  was  maximum.  (The  signifi¬ 
cance  of  this  quantity,  as  well  as  that  of  other  information  in  the  table,  is 
1  i  i- cussed  in  the  following  section.)  In  addition,  a  representative  set  of 
.  x-'ilf  te  experimental  alpha  particle  yield  curves  is  shown  in  Figures  12a  - 
1  .  it  should  be  noted  that  the  yield  curves  are  all  plotted  on  a  semileg¬ 

al'  ''  ;o  scale;  this  was  done  to  facilitate  visual  cross-comparison  as  well 
•  •  ophasize  the  importance  of  the  lower-lying  points  on  the  wings.  Super- 
'o  the  experimental  points  are  theoretical  vield  curves  which  were 
•  I  ;>c  ording  to  a  procedure  described  in  the  following  section. 


KEY  TO  FIGURES  12a,  b,  c  (pp  43,  44,  45) 

Observed  alpha  particle  yield  curves.  All  points 
shown  are  recorded  data  for  the  samples  which  are 
identified  at  left.  (Different  symbols  are  used 
for  ease  of  distinguishing  nearby  curves.)  How¬ 
ever,  for  the  lower  curves  (Samples  606,827,  A824 
and  802)  the  counting  time  was  doubled  to  obtain 
sufficient  statistical  accuracy.  The  points 
plotted  for  those  curves  are  half  the  recorded 
counts.  The  statistical  errors  (not  shown)  are 
the  square  roots  of  the  observed  counts;  for  the 
lower  curves  they  should  be  reduced  accordingly. 
Also  shown  with  each  observed  curve  is  a  theore¬ 
tically  derived  curve  (the  continuous  curve  be¬ 
tween  the  data  points) ;  these  were  selected  on 
the  basis  of  best  fit,  as  described  in  the  text. 
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TABLE  2 .  SUMMARY  OF  THERMAL  TREATMENTS  AND  RESULTS 


THERMAL  TREATMENT  PROFILING  RESULT  ANALYSIS 


* 


Sample 

Temp 

Duration 

Atmos¬ 

Peak  c4 

Fraction 

n 

D 

II 

(°C) 

(s) 

phere 

Count 

of  Ref. 

(by  peak 

frac)  (cm^/s) 

304 

413 

71640 

air 

3800 

0.970 

1.105 

1.23(10)"17 

308 

627 

22500 

air 

13674 

0.863 

1.49 

-16 

2.1 

309 

627 

230805 

air 

6491 

0.410 

4.85 

3.90u,/"16 

311 

788 

68521 

air 

2992 

0.473 

3.88 

8. 20 (10 )~ 16 

313 

789 

12600 

air 

5068 

0.801 

1.73 

6.27(10)"16 

316 

706 

30630 

air 

5144 

0.665 

2.32 

5.72(10)"16 

318 

707 

3540 

air 

7304 

0.965 

1.120 

2.87(10)'16 

401 

563 

1799210 

N2 

4882 

0.822 

1.65 

3. 79  (10)~18 

402 

548 

619400 

N2 

5155 

0.834 

1.58 

9.66(10)"18 

404 

543 

244900 

N2 

5634 

0.917 

1.295 

1.1K10)-17 

606 

609 

8520 

steam 

3005 

0.250 

8.10 

3.03(10)~14 

612 

609 

3111 

steam 

3111 

0.541 

3.12 

1.88(10)~14 

627 

1160 

11400 

N2 

4346 

0.785 

1.79 

7 . 73 (10)~16 

031 

690 

230700 

HV 

10909 

1 

1 

0 

631 

1150 

10980 

N2 

4353 

0.812 

1.68 

6.64(10)~16 

632 

815 

582240 

HV 

5087 

1 

1 

0 

707 

805 

7200 

Ho0 

3473 

0.721 

2.07 

1.82(10)~15 

716 

828 

(8.9  torr) 
253800  H2O 

10805 

0.719 

2.08 

5.27(10)~17 

717 

740 

(0. 

75900 

24  torr) 

h2o 

12781 

0.851 

1.54 

7.22(10)~17 

v  i  a 

81  8 

(0. 

16800 

80  torr) 

h2o 

8244 

0.865 

1.49 

2 . 87 (10)~16 

(2.65  torr) 
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TABLE  2.  SUMMARY  OF  THERMAL  TREATMENTS  AND  RESULTS  (contd.) 


THERMAL  TREATMENT  PROFILING  RESULT  ANALYSIS 


* 


Sample 

Temp 

Duration 

Atmos¬ 

Peak 

Fraction 

n 

a 

(°C) 

(s) 

phere 

Count 

of  Ref. 

(by  peak  frac)  (cm  / s) 

722 

842 

609000 

H20 

7621 

0.800 

1.73 

1. 30 (10)-17 

726 

514 

2700 

(0.062  torr) 
steam  7593 

0.812 

1.68 

2.70(10)'15 

/  28 

423 

5280 

steam 

7498 

0.786 

1.79 

1.67(10)'15 

7  3  .1 

331 

23640 

steam 

17589 

0.940 

1.213 

7.80 (10)-17 

263 

85440 

steam 

7573 

0.794 

1.75 

9.66(10)-17 

•<02 

476 

22380 

steam 

5828 

0.410 

4.80 

3.94(10)"15 

806 

476 

2775 

steam 

8899 

0.817 

1.66 

2.53(10)”15 

;-;n  i 

476 

7420 

steam 

9858 

0.706 

2.14 

1.92(10)"15 

All 

345 

25200 

steam 

11860 

0.850 

1.61 

2 . 53  (10)-16 

AS  15 

473 

7620 

steam 

12175 

0.901 

1.35 

4.36(10)-16 

816 

473 

7620 

steam 

9108 

0.674 

2.28 

2.20(10)"15 

A820 

475 

22770 

steam 

4905 

0.743 

1.97 

5.06(10)"16 

823 

619 

1290 

steam 

7919 

0.569 

2.90 

2.30(10)"14 

>.•24 

616 

16260 

steam 

5844 

0.455 

4.12 

3.93(10)-15 

.  \  >  L.  0 

619 

1290 

steam 

11246 

0.817 

1.66 

5.44(10)-15 

827 

616 

16260 

steam 

2419 

0.188 

11.80 

3.40(10)-14 

K:  Prefix  "A"  signifies  prior  high-vacuum  (HV)  anned:  1  hr  @  900  C. 
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Figure  12b.  Observed  alpha  particle  yield  curves  for  annealed  samples. 


ALPHA  PARTICLE  COUNT  PER  50uC  OF  INCIDENT  PROTONS 


.  i\\  ANALYSIS 


t  hoc!  of  Determining  Tracer  Diffusivitv 


Tin.1  tracer  diffusivity  D*  was  extracted  from  the  data  in  the  context  of 
.  w..:er  diffusion  model  discussed  in  the  Introduction.  According  to  this 
lit.'!,  the  exchange  of  network  oxygen  with  molecular  water  dissolved  in  the 
is  solely  responsible  for  anv  observed  change  in  the  distribution  of  im- 
i  ■'■!  ^  0.  The  dependence  of  D*  on  thermal  treatment  parameters  is  a  diag- 
:  of  the  mechanism  of  water  diffusion  in  Si O2  lavers  and  a  test  of  the 

|  1  i t v  of  this  model.  The  extent  of  change  in  the  distribution  is  in  prin- 
•  terminable  by  inversion  of  the  integral  eauation  (!4)  which  gives  the¬ 

ft  i c '■  e  vield  observed  by  nuclear  resonance  profil  ing.  In  practice, 

'  i or:  is  a  complicated  procedure  entailing  large  errors  whose  propa- 
1  : 1  u 1  *  to  t  rack . 

1 8 

"••ever,  once  tile  initial  0  distribution  is  determined,  its  subseouent 
is  predictable  through  macroscoo ic  diffusion  theory,  provided  the 
os  of  the  medium  are  correctlv  perceived  and  the  appropriate  boundary 
-s  are  applied.  A  full  mathematical  treatment  of  the  macroscopic  dif- 
•  1  awcss  governing  the  evolution  of  an  implanted  '  distribution  ap- 
Atoendix  A.  Tt  is  shown  there  that  when  a  sample  is  maintained  for 
at  a  fixed  temperature,  the  resulting  distribution  is  a  fun  r- 

single  parameter  -  namely,  the  product  D*t .  The  characteristics 
'‘fusion  process  are  embodied  in  this  parameter.  For  example,  its 
!»  ■  and  pressure  dependence  reflect  the  activation  energy  and  identi- 
: he  transported  species.  Its  position  or  time  dependence  indicate  non- 
ilies  in  the  medium  (see  the  Discussion).  For  the  case  of  diffusion, 
v,  the  problem  of  determining  0(%,  t)  is  reallv  that  of  evaluating 


"]'!v,  rather  than  parametrizing  the  distribution  with  D*t,  which  has 
•  - : on  a  dimensionless  parameter  was  chosen  for  claritv  of  presen¬ 

ts  parameter,  denoted  O,  can  he  thought  of  as  a  natural  time  unit 
: '  "e  the  extent  of  diffusion  in  terms  of  the  amount  of  broadening 
:>v  an  initial  Gaussian  distribution.  Its  applicabll  itv  to  this 
‘  ror  the  verv  close  similarity  of  the  implanted  Ho  concentration 
a  Ic’.miss  i  an  distribution.  it  is  well  known  that  a  Gaussian  distri- 
ivii'g  an  initial  variance  *i  broadens  with  time  in  an  infinite  homo- 
!  according  to  the  relation 

;’D*t  =  o?(t)  -  a." 

1 


is  variance  or  the  distribution  at  time  t.  t  c  H  be  defined 

.  tt  io  ,  i.e.,  the  extent  of  profile  broadening  of  a  Gaussian  dis- 

ion  tor  the  particular  case  of  an  infinite  homogeneous  medium.  However, 
can  he  used  to  index  the  broadening  of  a  Gaussian  implant  for  anv  other 
uni.irv  conditions-  The  corresponding  value  of  D*t  is  then  expressed  in  terms 
-  n  accord  ing  to  ^  =  (  2  _  })  ? 

1 


(HI 


4fS 


provided  I)*  is  independent  or  %  ana  t.  i'or  other  boundary  conditions  n  ioses 
its  meaning  as  the  ratio  of  widths,  but  still  serves  as  a  single  index  for 
the  product  D*t. 

To  determine  the  proper  value  of  c  to  be  associated  with,  a  particular 
thermal  treatment,  a  set  of  possible  discrete  values  for  n  was  selected.  For 
ea.  h  member  of  this  set,  a  distinct  distribution  p(X,t)  was  calculated  using 
the  method  described  in  Appendix  A.  Based  on  this  distribution,  a  predicted 
alpha  particle  yield  curve  Y(Ep,)  was  calculated  using,  the  method  described 
in  Appendix  B.  The  resulting  familv  of  yield  curves,  parametrized  bv  n.  was 
compared  to  the  experimental  vield  curve  resulting  from  that  thermal  treat¬ 
ment.  The  appropriate  value  of  n  was  taken  as  that  which  indexed  the  calcu¬ 
lated  yield  curve  most  closely  matching  the  experimental  vield  curve.  This 
procedure  allowed  the  extraction  of  the  relevant  information  on  diffusion 
without  resorting  to  inversion  of  the  integral  equation.  A  representative 
family  of  evolved.  l°0  concentration  profiles  calculated  in  this  manner  is 
shown  in  Figure  11,  while  a  familv  of  calculated  alpha  particle  vield  curves 
corresponding  to  these  distributions  is  illustrated  in  Figure  14, 

It  si. (hi.’  be  mentioned  that  a  number  of  preliminary  computations  were 
involved  i-.  th.  .-..1  euluc  ion  of  the  vield  curves.  These  involved  determining 
the  correct,  values  of  the  constants  appearing  in  the  formulas  for  the  reac- 
t ion  cross  section,  the  nenn  energy  sp  read  and  t  he  initial  J,0  distribution 
i  :i  Appena  :  A  .nd  ft.  A  trough  approx  in  .  te  v.i.uos  of  the  constants  were 

avail  an  t  e  r.n.  Mie  sour,  cs  cited  there,  tin-  finai  values  were  determined  in’ 
finding  t  rie  best  natch  of  the  calculated  initial  vield  curve  (topmost  curve 
in  Figure  id)  to  the  observed  vield  curve  of  an  implanted  but  otherwise  un¬ 
treated  simple.  The  goodness  of  fit  can  be  judged  bv  inspection  of  the  top¬ 
most  curve  in  Figure  12a. 

in  actual  practice,  the  procedure  for  determining  the  value  of  f\  rorres- 
.vndin ..  o  ,  .liven  t  ho  ran ;  treatment  consisted  of  taking  an  experimental 
vie!.;  .  ,  comparing  it  to  the  familv  of  precalculated  vield  curves,  then 

my  or  ram  vita  n  adjusted  to  obtain  as  good  a  match  as  possible. 

A .  •  iiou,..'.  i  not  niffirult  to  determine  n  with  an  accuracy  sufficient  for 

a  .  ,  -at  e  ,  nor gv  estimation,  tie  procedure  proved  unwieldv  when  higher  ac- 
.u/ecie,.  *"i  -  desired,  especially  for  small  values  of  n,  when  the  relative 
>  rr»,r  ir.  ..■  ■  trising  from  errors  in  n  grew  verv  large.  As  car.  be  seen  bv 
different .utlng  equation  (13),  the  relative  error  in  D*  is 

An*  _  2rT  An.  ( 16) 

p.f  n'l  _i  n 

For  large  values  of  n,  the  relative  error  in  D*  is  twice  that  in  •>.  As 
decreases,  however,  it  becomes  much  larger;  for  H  <  1.4  the  ratio  of  relative 
errors  is  -4.  In  addition,  the  shape  of  the  vield  curves  changed  relatively 
little  for  small  n,  making  the  determination  of  n  difficult. 

fo  overcome  these  problems,  an  alternate  procedure  was  adopted  for  rou¬ 
tine  me.  c  remen ts  of  it*.  This  procedure,  based  on  the  rate  of  decav  of  the 
peak  ainu..  particle  vield,  proved  especially  useful  in  the  regime  of  small  n. 
I:  .•an  •*.'  -,eer  f rer.i  the  familv  of  vield  curves  in  Figure  14  that  the  peak 
1  :  .  !  pli.i  , '.articles  is  a  monot  on  i  ea  1  1  v  decreasing  function  of  ft,  vThis 

function  is  shown  in  Figure  13.1  To  determine  the  value  of  n  using  this 
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14.  Theoretical  alpha  particle  yield  curves  (labels  indicate  n) 


Figure  15.  Time  behavior  of  theoretically  derived  peak  alpha  yield.  This  function  was  calculated  by 
calculating  a  large  number  of  theoretical  alpha  particle  yield  curves.  The  peak  value  of  each  of  those 
curves  was  tabulated  along  with  the  value  of  n  associated  with  the  curve.  The  figure  was  gotten  by 
connecting  those  values  with  a  smooth  curve.  The  variation  with  time  is  here  expressed  as  a  variation 
with  n,  as  per  the  discussion  in  the  text. 


function,  it  was  necessary  to  measure  an  alpha  particle  yield  curve  for  a 
sample  only  in  the  vicinity  of  its  peak,  so  as  to  determine  the  peak  yield. 
This  quantity  was  then  normalized  to  that  of  an  implanted  but  untreated  re¬ 
ference  sample  located  nearby  on  the  wafer.  The  appropriate  value  of  n  was 
then  found  simply  bv  locating  the  resulting  peak  yield  ratio  on  a  graph  of 
that  function  (Figure  15).  This  procedure  had  several  advantages.  First,  it 
greatly  shortened  the  time  spent  in  nuclear  profiling*  Second,  it  allowed 
shorter  thermal  treatments.  This  can  be  seen  from  Figures  14  and  15.  For 
relatively  small  0  (i.e.  short  treatments),  the  shape  of  the  yield  curve  is 
insensitive  to  changes  in  n,  so  that  comparison  of  shapes  is  a  poor  means  of 
determining  n.  However,  the  maximum  slope  of  the  peak  yield  ratio  occurs  in 
just  that  situation,  so  that  its  use  is  most  advantageous  there*  Third,  it 
allowed  errors  in  D*  to  be  estimated  directly. 

To  demonstrate  this,  consider  the  statistical  error  in  D*  obtained  by 
this  means  for  the  case  of  small  n*  Let  the  number  of  alpha  particles  counted 
be  denoted  N  and  N  for  the  treated  and  untreated  samples,  respectively .  The 
peak  yield  ratio,  denoted  R,  is  then 


The  statistical  error  in  R  is  gotten  from  the  statistical  errors  in  N  and 
N0  taken  in  quadrature,  namely 

(AR)2  =  (3R  AN)2  +  (SR  AN  )2  (17) 

3N  3N  ° 

o 

The  errors  in  those  quantities  are  just  their  square  roots.  If  these  are 
substituted  in  equation  (17),  we  find  that  AR  is  given  by 

“  -J¥o  ~ + «>  •  <18> 

According  to  Figure  15,  the  initial  functional  dependence  of  R  upon  n  is  ap¬ 
proximately 

R  -  ,  n  -  i 

R  1  4.1 

This  is  now  differentiated  to  get  Aft  in  terms  of  AR,  which  in  turn  is  given 
by  equation  (18).  By  substituting  the  resulting  expression  for  AH  in  equa- 
tion  (16),  the  relative  error  in  0*  is  found  to  be  AD*  _  8 . 2 f ~R(  1+R)  (19) 

n*  =  tfzrj  ^ 

With  N  =  5000,  AD*  is  about  20 7  for  11=1.5,  decreasing  to  10)'  for  11=2.5. 

°  D* 

The  calculation  of  error  propagation  can  be  easily  extended  to  include 
other  sources.  For  example,  lateral  nonuniformities  in  implanted  *^0  dose 
can  be  represented  bv  including  another  term  in  equation  (17).  (In  this 
study,  the  practice  of  using  nearbv  samples  reduced  the  magnitude  of  this 
error  to  below  the  level  of  statistical  uncertainty.)  Errors  in  temperature 
measurement  are  represented  bv  including  (in  quadrature)  the  additional  term 
Ha  AT  in  equation  (19). 
kT  T 

The  errors  in  reading  temperature  were  about  2  C.  For  the  range  of  tempera- 
f  romS  5%S to  257!°  ^  f°  ^  *  C  in  steanl^  thls  represented  a  contribution  ranging 
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''■  'ij'iTnture  and  Pressure  Dependence  of  Tracer  Diffusivity 

he  temperature  dependence  of  D*  was  determined  as  follows.  First,  D*  was  cal- 
. .  cd  from  the  observed  data  using  the  method  of  peak  yield  decay  described  in 
the  previous  section.  The  results  of  these  calculations  all  appear  in  the  last 
column  of  Table  II.  For  each  of  the  three  atmospheres  (i.e.  steam  at  1  atm,  room 
air  and  flowing  dry  nitrogen)  the  values  of  D*  were  fit  by  an  equation  of  the 
Arrhenius  type 

D*  ■  D*  e Cot) 

with  a  temperature-independent  pre-exponential  factor  D*.  The  values  of  D*  were 
taken  as  measured:  they  were  not  adjusted  for  the  temperature  variation  of  gas 
phase  concentration.  The  fit  was  carried  out  by  the  method  of  unweighted  linear 

■•as:  squares  using  -  and  log^D*  as  the  independent  and  dependent  variables, 

respectively.  The  results,  which  appear  in  Figure  16,  along  with  the  individual 
v.  of  D*,  are 


D*  =  2.7  (10) 

steam 


-L0+.40 


0*  .  =  1.3  (10) 

room  air 


-12+.34 


D*drv  N2  =  1.8  (10)~13 


-  ( 16.9  +  1.3  kcal  /rno  1 ) 

RT 

cm  /s 

-  (15.5  +  K4  kcal /mol) 


RT 

(15.7  keal/mol) 
RT 


o 

cm“/s 


on“/s 


(20) 

(21) 

(22) 


mcicated  errors  are  derived  solely  from  the  variances  of  the  fitting  param- 

(1),  and  take  no  account  of  the  sources  of  error  mentioned  in  the  previous 

;no.  The  standard  deviations  in  the  fits  for  log  D*  are  0.232  and  0.135  for 

■  and  room  air  respectively,  which  imply  relative  errors  of  171%  and  136%  for 

>■  dry  nitrogen  data  points  fell  fortuitously  close;  their  error  was  not 

ted.  The  three  fitted  lines  shown  in  Figure  16  are  almost  parallel;  at 

that  representing  D*  .  lies  a  factor  of  93  below  D*  ,  while 

&  room  air  steam 

lies  a  factor  of  8.1  below  that.  Also  shown  in  the  figure  are 

-  l l rogen 

,r  !)*  for  a  few  samples  that  had  been  vacuum  annealed  prior  to  thermal 

steam 

:  .1 T.e  n  t .  These  values  (the  square  symbols  in  the  figure)  lie  a  factor  of  5.0 
'  those  for  unannealed  samples.  The  reasons  for  this  will  be  discussed  shortly. 

hi.-  pressure  dependence  of  D*  was  determined  by  calculating  D*  in  a  similar 
>r  the  samples  heated  at  low  pressure  (#707,  716,  717,  718  and  722).  As 
been  heated  at  slightly  different  temperatures,  their  values  of  D*  were 
’  •  to  a  common  temperature  (of  820  C)  by  assuming  a  temperature  dependence 
i  seen  for  steam.  The  resulting  values  of  D*  (per  unit  pressure)  are  dis- 
■n  Figure  17.  Also  shown  is  D*steam  at  1  atm  which  was  calculated  using 
on  u-0).  The  indicated  errors  are  those  arising  from  counting  statistics, 
simple  (#632)  was  maintained  at  815  C  for  579600  s  at  a  pressure  of  (10)-(>  torr; 


W.  Mendenhall  and  R.  L  Scheaffer,  Mathematical  Statistics  and  Applications 
(Duxbury  1973)  pp.  391-393.  *  " 
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no  difference  could  be  discerned  between  the  yield  obtained  from  that  sample  and 
that  from  an  untreated  reference  sample:  at  that  pressure,  D*  was  below  the  limit 
of  measurement.  Another  sample  (#631)  had  its  alpha  particle  yield  measur  'wice. 
After  being  maintained  at  690  C  for  230400  s  at  (10)~“  torr ,  the  yield  was  npared 
to  that  of  a  reference  sample:  no  difference  could  be  detected.  It  was  then  main¬ 
tained  for  10980  s  at  1150  C  in  flowing  dry  nitrogen  at  1  atm  and  the  yield  was 
again  compared  to  that  of  the  reference  sample.  This  time  the  results  implied  a 
value  of  D*  close  ta  that  of  a  sample  (#627)  which  had  undergone  a  similar  treat¬ 
ment  in  dry  nitrogen  but  had  no  prior  vacuum  treatment. 

The  alternate  method  of  determining  D*,  comparing  the  shapes  of  observed  and 
calculated  alpha  particle  yield  curves,  is  illustrated  in  Figures  12a-c.  Accom¬ 
panying  each  observed  yield  curve  in  the  figures  is  a  calculated  curve,  a  member  of 
the  family  depicted  in  Figure  14.  The  basis  for  selection  was  visual  judgment 
rather  than  numerical  computation.  A  number  of  points  concerning  the  observations 
and  results  must  be  mentioned.  First,  the  alpha  particle  counts  indicated  in  the 
figure  are  either  the  actual  observations  or  (for  #606,  827,  A824  and  802)  hall  the 
recorded  counts.  The  corresponding  statistical  errors  in  the  counts  are  between 
1 %  and  2.5%;  as  they  can  be  calculated  directly,  they  did  not  seem  worth  indicating 
explicitly.  Second,  some  variation  in  overall  normalization  was  encountered.  This 
was  caused  by  lateral  nonuniformities  in  implanted  dose  and  by  slightly  difterent 
settings  of  detector  angle  on  different  runs.  The  maximum  variation  was  on  the 
order  of  10%. 

The  values  of  n  associated  with  each  of  the  theoretical  curves  in  Figures  12a- 
c,  as  well  as  the  resulting  estimates  for  D*,  are  listed  in  Table  III.  The  listed 
values  of  D*  for  the  samples  shown  in  Figure  12a  are  comparable  to  those  obtained 
by  the  method  of  peak,  yield  decay  listed  in  order  of  magnitude.  The  origin  of  this 
discrepancy  is  addressed  in  the  Discussion  section  which  follows. 


Figure  16.  (Next  page)  Observed  temperature  dependence  of 
180  tracer  diffusivity  at  1  atm.  The  solid  lines  represent 
linear  least  squares  fits  to  the  data  shown;  the  resulting 
pre-exponential  factors  and  activation  energies  are  indicated 
below. 


KEY  TO  FIGURE 


SYMBOL 

ATMOSPHERE 

Do* (cm2/s) 

Q  (kcal/mol) 

X 

Steam 

2.7  (10)~10 

16.9 

+ 

Room  Air 

1. 3 (10)~12 

15.5 

0 

Dry  N2 

1.8(10)-13 

15.7 

NOTE:  Q  represents  steam  measurements  of  samples  that 

had  previously  undergone  a  high-vacuum  anneal. 
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TABLE  3.  COMPARISON  OF  METHODS  FOR  OBTAINING  TRACER  DIFFUSIVITY 


YIELD  CURVE  SHAPE  PEAK  YIELD  RATIO 


Figure 

Sample 

n 

2  , 

D*  (cm  /s) 

n 

D*  (cm^/s) 

1 2a 

316 

1.75 

2.7(10)"]6 

2.32 

5.7(10)"16 

612 

2.5 

1.1(10)'U 

3.12 

-14 

1.9(10) 

606 

15 

5.5(10)"U 

8.1 

3.0(10)“14 

12b 

A826 

1.08 

5.2(10)"16 

1.66 

5.4  (10)-15 

A824 

1.5 

3.1(10)"16 

4.1 

3.9(10)"15 

12c 

821 

1.2 

2.4(10)_16 

1.69 

l.o(io)-15 

823 

1.5 

3.9(10)“15 

2.9 

2,3(10)"14 

802 

1.5 

2.2(10)"16 

4.8 

3.9  (10)~15 
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DISCUSSION 

Absence  of  Ion  Bombardment  Induced  Damage  in  SIO9  Network 


The  measurements  of  tracer  diffusivitv  performed  in  this  study  can  be 
compared  to  others  made  in  Si02  films  and  membranes  only  if  the  properties  of 
the  material  used  here  match  those  others.  In  view  of  the  well-known  ten¬ 
dency  of  ion  implantation  to  induce  lattice  damage  in  crystals,  it  may  be 
suspected  that  the  ^0  implantation  altered  the  structure  or  properties  of 
the  SiOy  network.  However,  a  good  deal  of  evidence  was  collected  during  the 
course  of  this  study  which  indicates  that  the  implantation  process  does  not 
produce  enough  network  damage  to  appreciably  alter  the  oxide's  properties,  at 
least  as  far  as  tracer  diffusion  is  concerned. 

First,  the  amount  of  ^80  implanted  was  small  enough  to  avoid  gross  mor¬ 
phological  changes  brought  about  by  its  incorporation,  such  as  blister  or 
bubble  formation.  (None  were  seen.)  The  additional  oxvgen  implanted  in  the 
network  amounted  to,  at  most,  1.9%  of  the  normal  concentration.  This  can  be 
seen  from  equation  (A4),  which  shows  the  peak  concentration  of  implanted 
in  the  initial  distribution  to  be 


r-  $ 

-J  2TT  - 


where  $  is  the  fluence  of  implanted  * 80o  ions.  As 
ions/cm^  and  was  0.0283  urn,  ^peak  was  4.2(10)20 

was  1.9%  of  the  concentration  of  0?  in  Si02»  which  is 
as  2.25(10)22  molecules/cm^ . 


this  was  3(10)1^ 
mol ecul es/rm 1 .  This 
given  in  equation  (7) 


Second,  the  shape  of  every  alpha  particle  yield  curve  could  be  well 
matched  to  one  or  another  of  the  set  of  theoretical  curves;  these  were  calcu¬ 
lated  on  the  basis  of  a  dilute  homogeneous  medium.  There  was,  therefore,  no 
indication  ever  observed  of  a  region  of  anomalous  diffusivitv. 

Third,  similar  thermal  treatments  were  performed  on  two  samples,  one  o' 
which  had  previously  undergone  a  high-vacuum  anneal  at  690  C  for  230400  s. 

The  180  tracer  diffusivity  was  found  to  be  unaffected  bv  this  anneal  .  These 
observations  can  be  interpreted  in  the  context  of  previous  invest igat ions  of 
annealing  of  radiation  damage  in  Si02  discussed  in  the  Introduction.  From 
these  studies  it  may  be  concluded  that  whatever  damage  is  produced  in  SiOj  is 
entirely  removed  by  a  suitable  high-temperature  anneal.  Tt  is  reasonable  to 
conclude  that  the  thermal  treatments  performed  to  induce  diffusion  them¬ 
selves  served  to  anneal  out  any  damage  induced  bv  the  implantation. 

Next,  there  is  the  question  of  whether  the  proton  bombardment  involved 
in  the  nuclear  resonance  Profiling  caused  anv  alteration  of  the  network  or 
induced  any  diffusion  of  l80.  To  answer  this,  alpha  particle  yield  curves 
were  collected  in  implanted  but  otherwise  untreated  reference  samples  both 
before  and  after  lengthy  proton  bombardments  (about  2000  uC,  at  a  current  of 
0.3  11A)  .  No  difference  could  be  detected. 
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Finally,  to  determine  whether  the  implanted  '  y0  actually  replaced  exis¬ 
ting  oxvgen  in  the  network,  or  remained  interstitial  am’  free  to  migrate,  an 
a'nha  part  ie  le  vie  Id  curve  was  taken  for  an  untreated  ref  erenee  samp !  e  attei 
a  b  I  gh- vacuum  anneal  at  XI  5  ('  for  S?9hnn  s.  jf  an  appreciable  amount  o'  im¬ 
plant  ed  *  s,t  had  remained  mobile,  then  it  would  have  d  i  '  •  used  during,  that  per¬ 
iod  with  <1  diffnsivitv  which  was  ciiarae  t  er  i  st  i  r  of  dissolved  molecular  oxveen 

■  — u  1 

at  that  temperature,  name  1  v  1  10)  cm'/s  (see  ,?ivurc  J).  !  lie  va  1  tie  o»  t  n> 

<.1  i  t  !  us  ion  iwr  iini'ter  n  corresponding  to  those  values  of  D*  and  t  can  he  c:,ci- 
l.tLed  using  count  ion  M  4  i  to  be  Iftno,  implying  the  escape  of  virtually  a!'. 

.'to  bile  0  through  the  oxide  surfaee.  However,  the  observed  yield  eurve  was 
i  ’nl  i  st  i  tutu  i  shah  I  e  '  *'>m  that  of  an  untreated  unannealed  referettee  sample,  i  u- 
d  (ratine  no  shift  in  the  initial  concent  rat  ion  profi  j.  Both  curves 

r  oselv  matched;  the  t  iieoret  ieal  1  v  predicted  curve  for  sacti  an  implant.  his 
deim  >nst  rat  e-:  t t  he  initia'  distribution  of  ^0  formed  !>v  the  imn !  an :  at  >  ■ 

••  '  sind  network  oxygen;  an  insignificant  fraet  i c.n  •• ; 

d:  nc !  ••ston  is  in  accord  with  the  stud  ie-  of  V 

o  •  :  •  v  mi  Si'  *•»  !>v  Hughes  (Johnson  et  a’- 33)}  .rui  ' 1  *  va  *•  •  •  •  .  ' 

•  '  <■  let  ro'.iuct  ion .  It  is  also  in  accord  with  :'p'\ *d*  • 

no.ivv  ion  energy  1  t,  which  nredids  bond  c  '  e  ■  va.'e  ■' 

. .  •  ■■■’■■ :  1  evolving  peak  tempera.  s  of  1150  d  :,f  •  pet  fpe  ,  - 

■  -  i  ■  .  ■  ■ :  ’  >•  .  c<-r  ion,  on  a  time  scale  of  (K‘  «  .  In  ' 

•• ,  *  -  -ax  i  mum  concentration  of  implanted  yn  Is  'esc  <  .■ 

'  . .  c.meentration.  The  subsequent  rat'id  ■'  •  •  - 

i  i  "  s-  . . '  etc  incorporat  ion  of  implanted  a-  ■  •  •  >u  «  •  '  . . he 

er  •  '’us  i on 

‘  -  discussed  in  the  previous  section,  the  implant  at  ion  •ausee  ’ v(!  to  re- 
a-  e  existing  oxvgen  as  a  constituent  or  the  XfO  ,  network.  'n  the  .d'sence 
water,  it  was  found  to  be  totally  immobile  at  temperatures  exceeding  H00C. 
vv--r,  in  *■  tie  presence  c>f  water  it  was  observed  to  move  in  a  manner  consi  s- 
'  h  !  t  s  'xehango  with  dissolved  molecular  water  an  1  its  subseeuent  d  i  r~ 
:•  is  a  (.  oust  it  lien  t  of  that  species.  'here  are  'out  piec<  s  o'  evidence 
•support  this  conclusion: 

The  observed  linear  dependence  of  the  0  tracer  di f fusivi tv  upon 
i"'  vap'-r  phase  concentration  is  identical  to  that  expect  ed  for  dissolved 
;  •  •  t  u'ar  water.  This  linear  dependence  was  seen  to  hold  over  4  orders  o'' 

.  »>;i  •  t  ude  (see  Figure  17),  fr  atm  to  79  ppm  (of  1  atm) .  These  measure- 
.  ats  were  made  in  pure  water  vapor  -  no  carrier  cases  were  employed.  Mien 
■  -  iir  wan  used  as  a  carrier  gas,  the  observed  tracer  diffnsivitv  conformed 

V  ■,  'inear  dependence,  but  depended  upon  the  partial  pressure  of  atmos- 
u;c.cr  vapor.  This  can  be  seen  from  Figure  lb  where  the  measured  ratio 

0*  steam  _  g ^ 

D*  room  air 

33.  N.  M.  Johnson,  W.  C.  Johnson  and  M.  A.  Lamport,  J  AppI  nhvs  46, 

1216  (1975). 

34.  D.  J.  DiiMaria,  D.  R.  Young,  W.  R.  Hunter  and  C.  M.  Serrano,  IBM  J  Res 
Devel ,  22,  289  (1978) 

93.  W.  Primak,  J.  Appl  Phvs  43,  2745  (1972) 
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is  close  Co  the  ratio  o?  water  v..:«>r  partial  pressures  for  these  atmospheres, 
namely  87  (92 
pies  were  tr«. 
after  6.7  davs  at  815C. 

(2)  The  temperature  dependence  of  D*  .-'■earn  showed  an  activation  energy 
of  16.9  +  1.3  kcal/mol,  which  is  indistinguishabl  >  ^rom  those  observed  for 
water  diffusion  bv  all  previous  investigators  (s.  igures  2,  3,  and  16). 
Previously  determined  activation  energies  range  from  16  to  18  kcal/mol.  This 
value  is  clearlv  inconsistent  with  those  previouslv  observed  for  oxvgen  dif- 
fusi^  ,  the  most  reliable  of  which  are  about  28  kcal/mol. 

(3)  The  magnitude  of  D*  is  close  to  that  expected  for  free  exchange  re¬ 

tween  network  oxvgen  and  dissolved  water.  In  the  presence  of  free  exchange, 
D*  is  related  to  the  experimentally  measured  diffusivity  of  water,  ac¬ 

cording  to  equation  (10)  bv 

/  d*  \  f  15  (10)”"  A00< 

V  Dfeff/  free  exchange  \  7.5(10)  “  '  ;  2(jflC  . 

This  ratio,  whose  value  is  derived  from  the  water  solui-i".  itv  measurements  of 
Moulson  and  Roberts  (54)  can  also  he  evaluated,  bv  using  the  results  of 
D*steam  obtained  in  this  studv  t see  Figure  Ini.  This  is 

D*steam  (this  work) 

Deff  (Moul son-Roberts) 

The  magnitude  of  D*<,team  is  thus  seen  to  he  consistent  with  free  exchange  at 
hoth  temperatures  measured  hy  Moulson  and  Roberts  to  within  a  factor  of  2-3, 
which  is  within  the  limits  of  accuracy  obtained  for  the  D*  measurements  in 
this  studv. 


J  (i  .  ii  \  :  0  )  t  h00( 

{4.4(10)““  •  1 20CH  . 


)  assuming  room  air  at  gOr  and  30?  rel  hum.  And  when  the  sam- 
'  1 8 

:ated  in  high  vacuum,  no  0  diffusion  could  be  detected  even 


(4)  In  the  presence  of  atmospheric  oxvgen,  the  activation  energv  for 
D*room  air  remained  that  of  water  (roughlv  16  kcal/mol,  as  indicated  in 
equation  21),  rather  than  revertinp  to  that  of  oxvgen  (28  kcal/mol).  This 
was  true  even  at  800  C,  where  exceeds  (Dpff)H20  bv  an  order  of  magnitude, 
and  even  though  the  concentration  of  atmospheric  oxvgen  was  17  times  that  of 
water  vapor  (20?  of  760  torr  vs.  50%  of  18  torr) .  This  implied  the  absence 
of  exchange  between  atmospheric  and  network  oxygen,  in  agreement  with  similar 
recent  findings  bv  Rosencher  et  al. (50).  As  discussed  in  the  Introduction, 
these  authors  observed  the  failure  of  diffusing  ^®02  to  replace  existing 
network  ^0  during  thermal  growth  of  oxide  films  on  silicon  in  dry  1^0  at¬ 
mospheres  . 


92.  R.  C.  Weast,  Handbook  of  Chemistry  and  Physics,  33th  ed  (CRC  Press,  1974) 
p.  D-159  '  .  -  ■  -  •  - 

54.  J.  Moulson  and  J.  P.  Roberts,  Trans  Faraday  Soc  5_7,  1208  (1961) 

50.  E.  Rosencher,  A.  Straboni,  S.  Rigo  and  H.  Amsel,  Appl  Phvs  Lett  34, 

254  (1979) 
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Water  I’ermeation  Kinetics 


The  values  of  D*  which  appear  in  Figure  16,  although  conforming  closely 
enough  to  an  Arrhenius  relation  to  support  the  conclusions  reached  in  the 
previous  section,  nonetheless  show  a  significant  mount  of  scatter.  This  is 
evident  from  the  large  standard  deviations  to  the  fits  o *'  equations  f 2 0 )  and 
which  are  6  times  larger  than  would  he  expected  on  the  >*asis  <«r  the 
discussion  of  equation  (19).  The  scatter  in  the  !>*  values  appears  to  ho  svs- 
t  em.tt  io  rather  than  random:  the  smallest  values  of  !»*  tend  to  he  as  sot  iated 
with  the  shortest  thermal  treatments  (i.e.  small  values  nt  n)  as  well  as  the 
!>••»  formaline  of  a  high-vacuum  anneal  prior  to  thermal  treatment  (so.  iahle  :  ' 
■  'irther,  the  amount  of  vield  curve  broadening  in  these  eases  is  .,n>  rta  1  •  ••.-.  " 
st.u  1  1  compared  to  the  decnv  in  peak  vield  (see  Figures  l.’  t-’J  i  .  these  -  ■  • 
'•at  ions  suggest  that  a  previous!  v  unsuspected  e* t  s  w; .  T',  t  •. 

v .  t  ’  "as  ,.f  !'*.  The  explan.it  ion  ‘or  the  s 'tread  >•  •*  v:  ■  •  ■  -  .  -  '<  •v.td 
•  ho  tinrontrol  V<!  cuantitv  or  roa'tec.  water  ■••  •  .  :  - 


an  initially  drv  oxide  ’aver  is  placed  in  n-  a*  '  h  — 

:•  s  water  vapor,  molecular  water  enters  the  ex  i  r.c  :  ,• 

'  ,c'  and  distributes  i t se 1  r  within  the  '  •> ve r  tl>at  i  •  s  ••  •  «  -t  v  •  ; 

•  •  '■!  '  v  attains  its  saturation  value.  At  t  he  out  set  ■■■'  this  studv,  i  t  was 

'  •'wed  that  equilihriur  between  dissolved  "i.'h"'u'ar  water  tnd  reacted  water 
the  net'  erk  is  also  ciuicklv  estah'ished,  so  that  tin  react  e.i  vote-  .  nm  o"- 

•  ion  so  rapidly  attains  its  saturation,  value.  However,  i  *'  t  h«  encroach 
-  eoui librium  is  not  rapid  and  the  oxide  1  aver  is  initially  drv,  there  w : 1 1 

a  de lav  observed  in  diffusion  until  the  reacted  water  which  is  neces- 

* r v  for  its  transport  has  rear hod  a  significant  level.  Tne  anomalous! v 
■•nail  values  observed  here  after  short  thermal  treatments  and  vacuum  desic¬ 
cation,  especially  for  the  lower  temperature  steam  treatments,  suggest  that 
these  conditions  occur. 

In  order  to  oroperlv  assess  the  effects  of  nonequi 1 ihrium  permeation 
•  it  ions  on  the  evolution  of  ^0  distributions,  a  thorough  theoretical  and 
imontal  studv  is  needed.  Although  a  complete  studv  is  bevond  the  s cove 
■  >  c  work,  an  introductory  development  appears  here  in  Appendix  D,  where 
■  riaatial  differential  equations  governing  1^0  diffusion  (a  formidable  set 
.*  v o  coupled  nonlinear  equations)  are  presented.  Thev  are  solved  here 
•'n ! v  ror  two  special  cases,  both  of  which  involve  rapid  saturation  of  the 
-xide  with  dissolved  molecular  water.  One  of  them,  involving  ran  id  estab¬ 
lishment  of  equilibrium,  is  the  case  which  was  thoroughly  treated  in  the  bodv 
•  this  study.  For  the  other  case,  the  content  of  reacted  water  remains  uni- 
..—i  throughout  the  oxide  but  varies  with  time.  Fcr  times  on  the  order  of  a 
•!  or  istir  time  T  it  is  far  lower  than  its  saturation  value,  causing  the 
'■•rved  initial  depression  in  D* . 

rt  is  interesting  that  this  same  phenomenon  seems  to  have  been  observed 
v  both  Will  lams (53)  and  Haul  and  Dumbgen  (51)  in  their  earlv  studies  o’  ''0 
permeation  in  SiOo  fibers  and  membranes.  As  discussed  in  the  Introduct ion , 

3.  R.  H.  Doremus,  Class  Science  (Wiley,  1971)  p.  12R 
53.  E.  Williams,  J  Am  Ceram  Soc  Aft,  190  (1965) 

51.  R.  Haul  and  0.  Dumbgen,  7.  Electrochem  66,  636  (1962) 
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these  authors  employed  mass  spectrometers  to  measure  atmospheric  oxygen  iso¬ 
topic  ratios  in  containers  with  limited  amounts  of  ^O-enriched  oxygen. 
Williams,  observing  that  the  values  of  diffusivitv  at  short  times  were  much 
lower  than  those  evaluated  at  long  times,  particularly  for  low-temperature 
diffusions,  invoked  this  effect  to  account  for  the  very  high  activation  ener¬ 
gies  found  earlier  by  Haul  and  Dumbgen  and  by  Sucov  (52),  whose  observations 
all  involved  short  times.  By  not  carrying  out  long  enough  diffusions,  they 
observed  erroneously  small  D  values  at  lower  temperatures,  leading  to  an 
overestimate  of  the  activation  energy.  The  activation  energy  he  found 
(29  kcal/mol)  was  close  to  that  measured  by  Norton  (49)  as  well  as  to  the 
thermal  oxidation  B  values  for  drv  oxygen.  The  magnitudes  of  the  diffusivitv 
measured  in  the  early  permeation  experiments,  however,  were  on  the  order 

of  those  measured  in  this  study  rather  than  those  of  Norton  and  the  thermal 
oxidation  experiments.  This  suggests  that  thev  were  observing  some  sort  of 
exchange  of  diffusing  ^0  with  network  oxygen.  Both  Meek  (941  and  Schaeffer 
and  Muehlenbachs  (95)  explained  the  six  order  of  magnitude  difference  on  the 
basis  of  a  direct  exchange.  However,  in  view  of  the  lack  of  exchange  ob¬ 
served  both  bv  Rosencher  et  aL  (50)  and  in  this  studv  for  ^0  diffusion  in 
dry  Si02  films,  their  explanations  must  lie  doubted.  If  any  exchange  occurs, 
it  is  more  likelv  mediated  by  reacted  water  in  the  oxide.  It  is  certainlv 
the  case  that  significant  quantities  of  water  were  present  in  Williams'  ex¬ 
periment,  both  in  the  oxides  and  in  the  atmosphere.  He  studied  its  influence 
on  his  diffusivitv  by  vacuum  desiccation  of  the  SiCH,  and  found  that  this 

procedure  hot;;  increase*:  the  initial  Jelav  and  lowered  the  limiting  value  of 
the  diffusivitv.  The  influence  of  water  was  apparent  Iv  not  appreciated  '->£ 
the  other  earlv  investigators,  who  seem  not  to  have  taken  pains  either  to 
eliminate  or  control  it.  The  actual  mechanism  of  exchange  (if  it  exists'  be¬ 
tween  diffusing  molecular  oxvger.  and  reacted  water  remains  obscure;  and  the 
process  has  not  been  observed  in  thermal  oxidation  experiments  which  have 
involved  oxygen  atmospheres  containing  controlled  fractions  of  water  vapor . 
The  parabolic  coef  f  ic  ion;  s  measured,  both  bv  Deal  et  al  (67)  and  bv  Ota  and 
Butler  (71s)  are  consistent  with  independent  diffusion  of  water  and  exveer 
over  a  range  of  utr.osnher  i  water  concentrations  from  to  1007  bv  volume. 
This  conflicts  with  t  h.  •;  -creation  of  Irene  and  Ghez  (40),  that  atmospheric 
water  concentrat  i, ..  .  .v  as  IB  ppm  (in02)  markedlv  increased  the  narahol  ic 
coe  f f  ic ient . 
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The  analvsis  of  these  authors  vielded  systematically  different  values  of  B 
from  those  of  other  investigators  (see  Figure  3,  particularly  their  values  of 
By-)  at  780  C  and  of  B^O  at  893  C).  In  addition,  the  dependence  of  their 
parabolic  coefficients  upon  the  crystallographic  orientation  of  the  Si  sub¬ 
strate  suggests  a  nonstandard  definition  of  B,  which  is  normal Iv  taken  to  in¬ 
volve  pure  bulk  diffusion  rather  than  interface  effects.  It  may  be  noted 
that  although  Deal  etaL(67)also  report  a  large  increase  in  B  with  the  addi¬ 
tion  of  trace  amounts  of  H2~>  based  on  data  plotted  in  their  Figure  3,  a  re- 
p lot  ting  of  that  figure  with  (Bohserved  -  B02)  replacing  B0bserved  as  ordinate 
(plotted  against  volume  %  H2O  in  O2)  shows  this  quantity  to  depend  strictly 
I  inearl v  upon  water  concentration,  as  was  found  in  this  studv.  Figure  ]8 
here  shows  their  data  plotted  in  such  a  manner,  with  the  ordinate  shown  as 
percent  of  R  in  pure  water  vapor  in  order  to  compensate  for  the  temperat uiv 
h-nendence  of  Bh?0-  The  linear  dependence  is  evident. 

Another  example  of  the  strong  influence  that  water  exerts  on  diffusion 
°i0'>  occurs  when  silicon  wafers  are  thermally  oxidized  in  resistance- 
•  i  .«'•.!  furnaces.  The  evolution  of  water  from  the  walls  e‘*  the  silica  d  i  rru- 

■  ’ v:  tubes  caused  a  great  deal  of  variation  in  early  results  on  thermal  oxi¬ 
dation  kinetics,  as  was  noted  by  Revesz  and  Evans  (62).  These  authors  e1 ir- 
: ted  the  wall  contribution  by  rf  heating  their  samples  in  a  cool-walled 

■11, ice.  Subsequent  investigators,  realizing  the  significance  of  this  effect. 

■  ■Moved  various  means  to  inhibit  water  vapor  evolution.  For  example,  Mavo 
”’d.  Evans  (  5)  used  polvsilicon  diffusion  tubes:  the  film  of  SiO?  which 

•  •  rmet!  on  the  inner  wall  acted  to  getter  anv  water  vapor  present  in  the  at- 

■aphere.  The  extent  to  which  silica  furnace  walls  contribute  to  atmospheric 
"or  vapor  content  can  be  seen  from  the  work  of  Vakavama  and  Collins  (691, 
m  measured  oxidation  rates  for  silicon  in  dry  argon  corresponding  to  water 
■•pur  pressures  of  0.22  torr  at  850  C  and  0.35  torr  at  1000  C.  The  same  ef- 
'  t  can  he  seen  from  the  values  of  D*drv  ^  observed  in  this  study.  These 
W'-re  consistently  a  factor  of  onlv  8  below  those  of  as  mentioned 

•*  tr Data  Analysis,  corresponding  to  water  vapor  pressures  of  about  1  torr. 

•  ■‘■kstreaming  of  room  air  through  the  open  exhaust  end  of  the  diffusion  tube 

also  have  contributed  to  the  large  value  of  P*drv  ^  •  As  no  quantitative 

■  urenres  were  drawn  from  that  data,  this  possibility  was  not  investigated 
'  1  he r  . 
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A.  0.  Revesz  and  R.  J.  Evans,  d  Phvs  Chem  Solids,  30 .  551  (1969) 
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Figure  18.  Increment  in  the  parabolic  rate  constant  upon  addition  of  H20  to 
oxygen  atmosphere,  as  percent  of  B^2o  in  pure  water  vapor.  All  data  points 
are  from  Deal  et  al.  (67) >  values  of  Bh20  in  pure  water  vapor  (used  for  nor¬ 
malization)  are  from  Deal  and  Grove  (56),  adjusted  to  100  vol  %  water  vapor. 


67.  B.  F..  Deal,  D.  W.  Hess,  J.  D.  Plummer  and  C.  P.  Ho,  J  Electrochem  Soc 
125,  339  (1978) 

36.  B.  E.  Deal  and  A.  S.  Grove,  J  Appl  Phys  36,  3770  (1965) 


An  independent  confirmation  of  the  time  variation  in  tracer  dif- 

fusivitv  resulting  from  nonequilibrium  permeation  conditions  has  recently 
been  provided  by  Rigo  et  al. (96) .  Using  procedures  similar  to  that  of 
Rosencher  et  al. (50),  they  measured  the  permeation  of  180 

into  Si02  films  upon 

thermal  treatment  in  atmospheres.  They  found  that  the  1&0  profiles 

were  consistent  with  a  time-varying  D*,  especially  for  temperatures  near 
those  employed  in  this  study. 


E.  Rosencher,  A.  Straboni,  S.  Rigo  and  0.  Amsel,  Appl  Phys  Lett  34_ 

254  (1979) 

S.  Rigo,  F.  Rochet  and  A.  Straboni,  "An  Study  of  the  Oxygen  Exchange 

in  S102  Films  During  Thermal  Treatment  in  Water  Vapor,"  Presentation, 
1980  International  Conference  on  the  Physics  of  MOS  Insulators, 

Raleigh,  NC  17-20  June  1980. 
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CONCLUSIONS 


» 


The  results  of  the  ^0  tracer  diffusion  measurements  performed  in  this 
study  clearly  demonstrate  that  water  diffusion  in  Si02  films  proceeds  by  the 
transport  of  molecular  water  and  that  it  is  accompanied  by  a  strong  reversi¬ 
ble  reaction  between  molecular  water  and  network  oxygen.  The  observed  linear 
dependence  of  tracer  diffusivity  upon  partial  pressure  of  atmospheric  water 
vapor,  for  pressures  ranging  from  1  atm  to  as  little  as  0.06  torr,  indicates 
that  network  oxygen  is  exchanged  with  molecular  water.  The  activation  energy 
of  tracer  diffusion  remained  about  16  kcal/mol,  characteristic  of  water  dif¬ 
fusion,  over  a  range  of  temperatures  from  1150  C  to  260  C.  Further,  the  mag¬ 
nitude  of  the  tracer  diffusivity  was  consistent  with  the  presence  of  free  ex¬ 
change  between  network  oxygen  and  diffusing  water.  The  extent  of  tracer  dif¬ 
fusion  proved  to  be  unaffected  by  the  presence  of  atmospheric  oxygen,  indi¬ 
cating  the  absence  of  exchange  with  that  species.  The  results  obtained  here 
agreed  with  corresponding  results  obtained  by  most  other  investigators  using 
different  means.  They  are  in  complete  agreement  with  recently  proposed  mo¬ 
dels  of  oxygen  and  water  diffusion  in  SiO,,  in  which  the  diffusion  Droceeds 
by  transport  of  the  molecular  species  through  interstices  in  the  network 
without  any  direct  reaction  between  the  molecular  species  and  the  Si02.  (The 
strong  exchange  reaction  accompanying  water  diffusion  is  mediated  by  immobile 
OH  groups  which  are  dissociation  products  of  the  molecular  water.) 

The  process  employed  for  introducing  tracer  into  the  network,  that 

of  ion  implantation,  was  shown  to  produce  a  region  within  the  oxide  in  which 
the  ^0  was  fully  incorporated  in  the  network:  no^O  diffusion  was  observed 
in  high  va'cuum, even  after  161  hrs  at  815  C,  The  properties  of  the  network 
remained  otherwise  unaltered,  at  least  as  far  as  diffusion  was  concerned. 

This  was  demonstrated  both  bv  the  data's  close  match  to  theoretical  predic¬ 
tions  based  on  a  dilute  homogeneous  medium  as  well  as  by  the  absence  of  any 
effect  on  diffusion  of  a  high-temperature  treatment  which  was  sufficient  to 
anneal  out  anv  implantation- induced  damage.  These  findings  indicate  that  the 
thermal  treatments  performed  to  induce  ■‘■^0  diffusion  were  themselves  suffi¬ 
cient  to  anneal  out  any  damage. 

The  technique  of  nuclear  resonance  profiling  proved  to  be  a  convenient, 
nondestructive  and  sensitive  method  of  measuring  ^"0  diffusivities  as  low  as 
1(10)-17  cm^/s  at  temperatures  as  low  as  260  C. 

In  addition  to  supporting  the  above  conclusions  concerning  the  mecha¬ 
nical  water  diffusion  in  Si02,  the  measurements  made  in  this  study  provided 
evidence  indicating  that  chemical  equilibrium  between  water  dissolved  in 
Si02  films  and  its  reaction  products  often  had  not  been  attained  over  the 
time  spans  of  the  thermal  treatments.  The  time  variation  of  the  reacted 

1  O 

water  concentration,  leading  to  a  time  variation  in  iO0  tracer  diffusivity, 
suggests  a  series  of  ^0  diffusion  measurements  which  would  form  a  natural 
follow-up  to  this  study.  Thermal  oxide  films  could  be  grown  on  a  set  of 
silicon  samples  in  dry  oxygen,  each  of  which  would  contain  a  thin  layer  en¬ 
riched  with  (formed  either  by  a  short  oxidation  in^O  or  by  a  low-energy 

implant) .  Although  the  total  thicknesses  and  the  enriched-laver  thicknesses 
would  all  be  identical,  the  depth  of  the  enriched  layer  would  be  different 
for  different  samples.  After  the  set  of  samples  received  a  thermal  treatment 
in  water  vapor,  the  concentration  profile  of  reacted  water  in  the  oxide  could 
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tie  determined  by  measurement  of  D*  near  each  enriched  layer  by  nuclear  re¬ 
sonance  profiling.  Successive  concentration  profiles  could  be  determined 
from  a  series  of  thermal  treatments,  providing  a  complete  experimental  deter¬ 
mination  of  nonequilibrium  water  permeation  kinetics. 

Such  a  study  could  be  of  great  benefit  in  the  manufacture  of  future 
high-reliability  integrated  circuit  devices.  The  absorption  of  moisture  is 
known  to  seriously  degrade  the  breakdown  resistance  of  MOS  gate  oxides, 
leading  to  premature  device  failure.  Detailed  observations  of  permeation 
kinetics  would  be  a  valuable  tool  in  the  investigation  of  new  methods  to 
inhibit  water  permeation  in  oxide  films. 
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APPENDIX  A  -  Evolution  of  0  Distribution 

18 

Suppose  an  Si(>2  layer  containing  implanted  0  is  heated  for  a  time  t  under 
conditions  which  produce  a  tracer  diffusivity  D*  which  is  independent  of  x  and  t, 
where  x  is  depth  in  the  layer.  Then  the  distribution  of  excess  ^-^0  in  the  oxide 
(i.e.,  in  excess  of  the  natural  background  level  of  0_02%  isotopic  fraction)  will 
change  with  time  so  as  to  satisfy  the  diffusion  equation 

2  2 

/(x,t)  =  D*  /(x,t)  (Al) 

dx 

under  boundary  conditions  appropriate  for  the  experimental  case.  The  general 
solution  to  this  equation  is 

2  * 

=  A0  +  B0*+  5^~an  C  LAnC#5an*'+  BnSin  V**]  *  (A2) 


Remember  that  0  is  transported  only  as  a  constituent  of  molecularly  dissolved 
water.  The  boundary  conditions  are  therefore  determined  by  considering  the 
behavior  of  this  substance  at  the  boundaries: 

(a.)  At  the  oxide  surface,  equilibrium  between  the  outside  atmosphere  and  water 
dissolved  in  the  oxide  occurs  rapidly:  surface  absorption  and  reactions  are  not 
involved  (1).  This  implies  the  condition  ^(0,t)  =  0. 

(b)  At  the  oxide-silicon  interface,  the  boundary  condition  is  temperature 
dependent.  This  can  be  seen  from  the  analysis  of  steady-state  steam  oxidation 
performed  by  Deal  and  Grove  (2).  (Steady-state  means  V'  J  =  0  within  the  oxide) 
According  to  their  equations  (6)  and  (7),  the  ratio  of  molecularly  dissolved 
water  concentration  at  the  interface  to  that  at  the  outer  surface  is  given  by 


C  1  +  k*  7d  ~ 

o  o  ef  f 

where  k,  the  interface  rate  constant,  is  determined  through  observations  of  the 
linear  and  parabolic  rate  constants  A  and  B  and  the  equilibrium  solubility  C* 
(*o  is  the  layer  thickness.)  Deal  and  Grove  give  a  value  of 

k  =  1.8  (10)B  um/hr 

at  1000  C  with  an  activation  energy  of  45.3  kcal/mol.  Using  this  in  conduction 
with  Moulson  &  Roberts’  (3)  equation 

D  ff  =  (10)  B  exp  (-l8.3(kcal/mol)/RT)  (cm2/sec). 


A-l.  R.  H.  Doremus,  Glass  Science  (Wiley,  1973),  p.  128 

A-.?.  B.  E.  Deal  and  A.  S.  Grove,  J  Appl  Phys  36,  3770  (1965) 

A-3.  J.  Moulson  and  J.  P.  Roberts,  Trans  Faraday  Soc  57,  1208  (1961) 


ti-fc  temperature  dependence  of  is  seen  to  be  the  following: 

o 


T  (C) 

kX  ,D  ,, 
o/  eff 

h/co 

250 

3.1(10)~7 

i 

475 

7.7(10)"4 

i 

600 

0.011 

0.989 

1000 

1.4 

0.417 

1200 

5.9 

0.145  . 

temperatures  encountered 

in  this  work,  the 

r . 

i  . 

ratio  —  is 

'ate  boundary  condition 

is  therefore 

G 

o 

is  more  than  98/'.  The 


\  - / 
£  J- 

V-’  -• 


a-/)(*0,t)  .  0. 


18. 


lor  sufficiently  long  times  all  the  excess  0  will  eventually  disappear, 
through  surface  leakage  or  interface  incorporation  of  diffusing  molecular 
implying  the  condition  J>(x,o6)  =  Q. 


(c)  The  initial  distribution  was  taken  to  be  that  implantation  profile  given  by 
Gilbeiis,et  al  (4)  for  the  appropriate  experimental  conditions.  This  was  an  Edge- 
worth  distribution  (the  product  of  a  Gaussian  and  a  polynomial);  explicitly,  it 


_  .  Fluence 

/)(x,0)  =  -2---  g  2 


imp 


-iTi  +v4i  3  0  + 


/s 
6 


h 
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(C4  -  6C2  +  3)  +  "5-  “  15^4  +  45^2  -  15>1 

2  72  J 


(A4) 


x-x, 

E,  m  .  _JLSE 
0  imp 


X.  =  0,0768  jUm 
^imp 

a  =  0.0283  Aim 
imp 


/T. 


1  =  -  0.164  . 


A-l.  J.  F.  Gibbons,  W.  S.  Johnson  and  S.  W.  Mylroie,  Projected  Range  Statistics, 
2nd  Edn.  (Academic,  1975) 
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Conditions  (a),  (b) ,  and  (c)  imply 


A  =  B  =  A  =0  (all  n) 
o  o  n 

B  =  0  (n  even) 
n 


The  solution  to  the  diffusion  equation  is  then  straightforward,  being 

.  prir  ¥) 2  D*tl  .  /2u-i ,  \ 
,-ie  ■>  —  v 


-P(*»t>  =  £ 


n=i 


2n- 


with 


\ =  2^dx  sin  (¥  x)/(x’0>  • 


(A5) 


(A6) 


The  numerical  evaluation  of  J9(X,  t)  was  performed  with  an  Interdata  832  Computer 
located  at  Fort  Monmouth,  N.J.  The  programs  were  straightforward.  Representative 
results  appear  in  Figure  13. 
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APPENDIX  B  -  Evolution  of  Alpha  Particle  Yield 


18  IS 

The  yield  of  alpha  particles  from  the  0(p,  a  )  N  reaction  is  Riven  by 


the  integral  equation 


/*xo  r a’  E  • 

Y(Eb>  Vo  dx/(x’t)/odEi  8(Eb’Ei/0  1dE^r(E)f(E,Eitx) 


x )  ( B  I  ) 


where 

Eb  is  the  mean  energy  of  the  incident  proton  beam, 

E  is  the  energy  of  reacting  protons, 

E^  is  the  actual  energy  of  protons  incident  on  the  surface  of 

the  layer, 

x0  is  the  thickness  of  the  oxide  layer, 

is  the  depth  at  which  the  reaction  occurs, 

18 

is  the  0  concentration  profile  (i.e.  atoms  per  unit  volume 
at  depth  x  at  time  t)  , 

is  the  differential  cross  section  of  the  ^0(p,  ot  )E',\' 
reaction, 

g(EjJ,E^)dE^  is  the  fraction  of  protons  in  the  beam  which  have  an  incident 
energy  between  E^  and  E^  +  dE^, 

f(E,E..x)dE  is  the  probability  that  a  proton  which  has  started  with 

energy  E^  and  arrived  at  depth  x  will  have  an  energy  between 
E  and  E  +  dE. 

This  also  appears  in  the  Introduction  as  equation  (14).  Once  the  functions 

J 

/W),  g(Eb,Ei),  f (E.Ei.x)  and^y(E)  are  established,  then  the  yield  is 
obtained  by  integration.  A  means  of  calculating^/5 is  given  in  Appendix  A. 
For  the  remaining  functions,  we  proceed  as  follows: 

g(Eb,E1) : 


The  distribution  of  energies  in  the  incident  beam  was  found  experimen¬ 
tally  through  analysis  of  the  thick-target  Al(p,^)2®Si  excitation  curve 
(see  Procedure).  As  the  energy  spread  amounted  to  less  than  1  keV,  the 
details  of  its  shape  are  relatively  unimportant;  for  ease  in  calculation  it 
can  be  taken  to  be  a  Gaussian  with  °b&l  keV,  i.e., 

-(Eb-Ei)2 

8(Eb,Ei)  =  — —  e  - “  052) 

/2^ab  20  b 


X 

yO(x,t) 

d  o 
da 
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t <K,E±,x) : 


If  the  proton  energy  loss  were  strictly  deterministic,  then  f  would  be  given 

S/'X  ,dE  ,  dE 

(E-fE^-w^dx  —,,}),  where  — f  is  the  specific  energy  loss  of  protons  in 
X  X  dE 


SiO 


2_  At  an  E^  of  629  keV  the  value  of  is  61  keV/um  (1) 


Since  the 


layer  is  only  0.2  um  thick,  the  protons  lose  12.2  keV  in  a  complete  traver- 

dE 

sal;  throughout  this  region  -g^>is  constant  to  within  1%.  Thus 


,dE 
dx  — > 
o  dx 


0 


x  = 


ex 


However,  the  random  nature  of  the  independent  scattering  events  by  which 
protons  lose  energy  in  materials  implies  that  their  energy  at  a  given 
projected  range  is  not  uniquely  determined  by  their  initial  energy.  The 
effect  of  this  energy  straggling  is  to  broaden  the  delta  function  to  a 
distribution  f(x,A),  where  x  is  the  projected  pathlength  and  A=  E^-E(x). 
Vavilov  (2)  found  that  this  function  may  be  closely  approximated  by  one  of 
three  expressions,  the  choice  depending  on  the  magnitude  of  x.  (A  review 
of  Vavilov's  work  appears  in  Appendix  C„)  Upon  performing  detailed  calcula1 
tions  based  on  that  work,  it  was  found  that  for  the  particular  conditions 
of  this  experiment,  the  distribution  could  be  approximated  by 


f(x,£) 


1 

/2tt  kgX1 


e  *p 


where  kg 


17.75  (keV2  um-1) 


(A-<x)  2 
2ks  x 


(B3) 


dflr 

d« 


(E): 


Highly  accurate  measurements  of  the  differential  cross  section  for  the 
18  15 

0(p,a  )  N  reaction  at  Q  =  150°  have  recently  been  performed  by  Amsel, 
Maurel  and  Nadai  (3).  3 


B-L  E.  Bonderup  and  P.  Hvelplund,  Phys  Rev  A4,  562  (1971) 

8-2.  P.  V.  Vavilov,  Soviet  Phys  JETP  _5>  749  (1957) 

8-3.  J.  W,  Mayer  and  E.  Rimini,  Eds.,  Ion  Beam  Handbook  for  Materials  Analysis 
(Academic,  1977)  p  163 
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A  very  close  fit  to  their  data  in  the  region  of  the  629  keV  resonance  ha* 
been  found  to  be  a  sum  of  two  interfering  Lorentzians  plus  an  exponential 
background.  Explicitly,  this  was 


~  (12.78-K)  -  (E-Eres)  4.7£ 

(E-Eres)  +  [  1  p4 
1  res 


E-Eres 

116 


E-Eres 

+  K6  58  (mb/ste) 


where 


2  (E  <  Eres) 
1.75  (E  >  Eres) 


P  =  1.9  keV 
res 


E  =  629  keV  . 
res 


With  all  the  constituents  of  the  yield  established,  its  evaluation  was 
carried  out  numerically  by  first  substituting  equations  (B2)  and  (B3)  in 
equation  (Bl).  Because  of  the  strong  damping  of  the  exponential  function 
f(x,&)  in  the  range  of  %  where  P(x  , t)  was  appreciable,  the  limits  on  the  E 
integral  could  be  extended  to  -OO  without  .ignif icantlv  changing  the  result. 
Similarly,  the  strong  damping  of  the  exponential  function  g(Eb,Ep  allowed 
the  lower  limit  on  the  integral  to  be  extended  to  The  order  of 

integration  was  then  interchanged,  enabling  the  yield  to  be  written 


/°  ✓  s 130  r~  ”1-1 

Y(Eb>  =  /dx^(x,t)  I  dE  ~(E)  /  dEi  4n  Oh  v'k”  x  ' 


exP  n  - 


(Ei-Eb)t  (F-i-  [E+Ex]  ) 


2ks  x 


The  E^  integral  could  then  be  performed  using  the  convolution  identity 


f  O  (5-t) 

Ut  e‘»p  e' v 


0^  0,,  /tT* 

/  2  2~ 

fC,  +  0, 


0  C  P 
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with 


t  =  Ej_-Eb,  S  =  E-Eb+Cx,  Cj*  >/?  :b,C5=  /2ks 


The  yield  is 


Y(Eb) 


-  X 

dx 


o  00  (E-Eb+€x) 

°  >  /  dEt(E>e'^~rv^-  <b?> 

2 17  (°b  +  kSx>  /-c°  b  s 


,  E-(Eb  -  EX) 

Using  the  substitution  E  =/ - ^ -  ,  this  was  rewritten  as 


7^? 


+  k  x) 
s 


Y(E, 


,> -  ^  +  V>  +  E„  -«*)(b8) 


This  was  the  appropriate  form  for  evaluation  by  Hermite  integration  (4); 


/  dE*  f(E')  ^  W  f(E/,1)  . 

j— j  3 


The  evaluation  of  equation  (B8)  was  performed  using  the  weight  factors  and 

.  d_£ 

and  abscissas  appropriate  for  a  20-point  Hermite  integration  over  E  with 
given  by  equation  (B4),  and  a  50-point  Simpson  integration  over  x,  with 
>>(*’t)  given  by  the  sum  of  equation  (A5)  plus  a  constant  background).  A 
listing  of  the  program  appears  on  the  next  three  pages;  it  was  run  on  an 
Interdata  832  computer  located  at  Fort  Monmouth,  N. J .  Representative  results 
appear  in  Figure  13. 


B-4.M.  Abramowitz  and  I.  A.  Stegun, 
(US  Government  Printing  Office, 


Handbook  of  Mathematical  Functions , 
1964) 
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APPENDIX  C  -  Proton  Energy  Straggling  Distribution 


If  a  proton  with  initial  energy  E^  passes  through  a  medium,  its  energy 
E  decreases  by  an  amount  A  =  E^-E(x)  where  x  is  the  projected  pathlength. 
Because  of  the  random  nature  of  the  independent  scattering  events  by  which 
protons  lose  energy  in  materials,  the  proton  energy  at  a  given  projected 
pathlength  is  not  a  unique  function  of  E-j.  Rather,  the  proton  energy  loss  is 
characterized  by  a  distribution  f(x,A).  f(x,A)dA  is  the  probability  that 
a  proton  traversing  a  path  of  length  x  will  suffer  an  energy  loss  between  A 
and  A  +dA  .  This  distribution  was  studied  by  Vavilov  (1),  who  found  that  it 
is  closely  approximated  by  one  of  three  expressions,  the  choice  depending  on 
the  magnitude  of  x.  (The  equations  given  by  Vavilov  will  be  denoted  herein 
by  the  prefix  "V".) 


For  x  -  1  um,  the  expression  for  f(x,A)  given  in  equation  (V15)  is 
equivalent  to 


f(x’A)  =  v^TT7  exp 

S 


(A-fo)  ' 

2k  x 
s 


(Cl) 


wxtn 


,  2  ,  4 

ksx  =  “g  =  ^ ~  z^e  Nx 


where  fig  =  Bohr  straggling  parameter  (Bohr  (2)  as  quoted  by  Fano  (3)) 
z  =  proton  atomic  number 
Z  =  Si02  atomic  number 
e  =  charge  of  the  electron 
N  =  number  density  of  Si09  atoms 

a  <  < 

,'c  1000  A  -  x  -  1  um,  f(x,A)  is  given  in  equation  (V13)  as 

■  (x,  A  )  =  e^t_  f  )  Ai  ( t )  •  (C 2 ) 

Tiiis  expr  s  ion  can  be  evaluated  as  follows: 

Accord i ng  t  «>  the  ci event ional  notation  of  special  relativity,  the  protons' 
velocity  (re’  t i ve  to  the  speed  of  light)  is 


2.  2E 

*  =MC2 

C-l.  P.  V.  Vavilov,  Soviet  Phys  JETP  _5,  749  (1957) 

C-2.  N.  Bohr,  Mat.  Kys.  Medd.  Dan.  Vid.  Selsk  _1_8  (8)  (1948) 
r.-?..  !  .  l'ano,  Ann  Rev  Nuel  Sc  i .  1  3,  1  (1963) 


( C  3 ) 
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where  M  -  proton  mass  and  E  =  proton  energy.  From  equation  (V5)  we  have 

(C4 ) 


£  =  smalJ^^L  _4_ 

‘'max  ,  „  2  ..  2  1836 


1  -£■ 


Me" 


where  n„  =  electron  mass. 


And  also  from  equation  (V5) 

5  =  0,300  x 


™eC  Z 
A  2  A 


0.075  x 


meC2MC2 


(05) 


where  x  =  depth  within  layer,  in  g/cm,  and  z/A  =  1/2 

" rom  he- tween  (V5)  and  (V6) ,  and  using  equations  (C4)  and  (C5) ,  we  have 

„  ,2y.2 


K.  =  ~  =  34.43  x 


(C6) 


,2 


equation  (V12),  and  using  equation  (C6) 


o.-f2) 


2  K 


(1-2/3  (32)2 


1/3  [2*]1/3  =  4.099 


2  2 
xm  C  M  C~ 
e 


1/3 


( C  7 ) 


Also  from  equation  (V12),  and  using  equations  (C5)  and  tC7) 


V 


r  ■>  _ 

1/3 

(1-  3  P2) 

small  0  ^ 

-  =  0.004464 

r*  ~ 

2  2 

xm  C  MC-  E 
e 

a 

—  - 

1/3 


<C8) 


Expressing  x  in  uni,  with  p=  2.27  g/m3  for  SiOo  (4)  and  approximating  E  by 
Eros  =  629  keV,  we  have 


\  = 

0.01298  x 

(MeV) 

(C9) 

/C« 

9.472  x 

(d imensionl ess) 

(CIO) 

a  = 

2.666  xl/3 

(d  imens i on  1  ess) 

(Cl  1  ) 

and 

n =  1 .826  x5' 3 

(keV)  . 

'"C 1  2) 

C-4.  A.  S.  Grove,  Physics  and  Technology  of  Semiconductor  Devices  (Wiley,  1967) 
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r 


Using  equations  (Cll)  and  (C12)  and  £  =  61  keV/um  (5)  in  equation  (V12) , 
we  have 

(A-gx) 


t  = 


+  a* 


(C13) 


The  distribution  given  in  equation  (V13)  can  then  be  written  as 
0.309 


0  309  / 

f  (x, A)  =  ~~T J  exp(1.46A  -  96.43x)  Ai |  — 

x  \  1. 


A-'61  +7.108x2/3 


826x 


1/3 


(C14) 


where  Ai  is  the  Airy  function  (6) . 


The  following  give  an  excellent  fit  to  all  values  of  Ai(x)  listed  in 
Table  10.11  of  Reference  6: 


r  -1/4 

Ixt  '  (0.3975cos4  +  0.4003  sin£) , 


(x<  -2.5) 


A)  (%)  =  { 


0.35503  (l  +  \  x 
0.25882  (x 


1  ..3  4  6 

+  -zr—r  X  + 


720 


28  9 

x 


362880 


_2_  *  10  7  80 

24  X  5040  X  3628800 


)- 


(-2.5<x<  1) 


1  i/^-c 

2  * 


0.00481  ^1 . 5  -  j)  2  +  0.01756(1.5  -  7)+  0.52703^, 

(1<  x)  ’  (C15) 


where  i  =  7!* 


2,..,  3/2 


For  x^lOOO  R,  f (x, A)  is  given  by  equation  (V16) : 


f (x,A) 


1  Ke*1+fiC)j  dy  e*fl  cos(yAi  +  #Cf2) 


(C16) 


C-5.  E.  Bonderup  and  P.  Hvelplund,  Phvs  Rev  A4,  562  (1971) 

C-6.  A.  Abramowitz  and  I.  A.  Stegun,  Handbook  of  Mathematical  Functions,  446 
(US  Government  Printing  Office,  1964) 
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where 


c  =  0 . 577  .... (Euler ' s  constant) 

=  fi2  [  In  y  “  Ci(y)]  -  cos  y  -  y  Sl(y) 

f2  =  y  lln  y  ~  +  sin  y  +j62Sl(y) 

=  ICa  +  k.  In  jc 


-  1  -  fi2  +  C 


In  K. 


and  where  Ci  and  Si  are  the  cosine  and  sine  integral  functions  respective 

Although  later  calculations  of  f(x,A)  exist  (see,  for  example  Bichsel 
and  Saxon  (7)  and  references  cited  therein),  recent  straggling  experiment 
involving  0-2  MeV  protons  in  thin  films  of  silicon  (8)  have  confirmed  tVj 
validity  of  Vavilov's  expressions. 


C-7.  H.  Bichsel  and  P.  Saxon,  Phys  Rev  All,  128ft  (1975) 
C-8.  J.  Baglin,  private  communication  (1979) 
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APPENDIX  D  -  Nonequilibrium  Permeation  Kinetics 

18 

To  fully  describe  the  evolution  of  the  tracer  0  distribution  in  the 
SiOo  films,  it  must  be  realized  that  can  be  a  constituent  of  either  dis¬ 
solved  molecular  water,  reacted  water  or  network  oxide.  Let  the  respective 
concentrations  of  the  various  species  be  denoted 


I  A(x.O 

f  h216o 1 

= 

[«218»J 

'  (x,  t) 

r  16  , 

[  Si  OHJ 

l/'Ax.t) 

= 

r  18  , 

[Si  OH] 

r c*,t) 

i 

s 

-  ,  .  16  . 

i  i:’1  •* 

!/>.<> 

c  18  16 

[  Si  0- 

(dissolved  water) 


(reacted  water) 


(network  oxide) 


18 

(The  stars  indicate  the  presence  of  tracer  0.)  We  may  assume  that 

a.  the  amount  of  ^0  is  everywhere  small  compared  to  the  amount  of  0, 
and  thus  neglect  any  terms  involving  r*r- 

b>  ^^(x.t)  =  y^SiO^j  independent  of  x  and  t,  i.e.  a  uniform  oxide. 

Under  these  conditions  the  equations  governing  the  tracer  distributions  are 


V/ 


v>/  -  -  <7/  -  ) 


(y<) 


dfr 

D-t 


7 ( -111  * 


(Dl) 


Vo 


o,  eg  -il m 

Z  K'v 


(D3) 


at 


v 


(D4) 
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05) 


'('lit'  first  tern  in  equations  (Dl)  and  (D4)  describes  the  transport  of  dis¬ 
solved  molecular  water;  the  quantity  is  the  flux  of  that  species.  The 

other  terms  in  the.  equations  describe  the  formation  and  disappearance  of  the 
various  species  through  the  processes  of  recombination  and  dissociation; 
terms  involving  products  of  concentrations  describe  binary  collisions.  The 
uiantity  T  can  be  thought  of  as  the  lifetime  of  a  dissolved  water  molecule, 
ualogously  to  the  charge  carrier  lifetime  which  appears  in  the  equations 
describing  electron  and  hole  diffusion  in  solids  (1).  it  should  be  n<>ted 
that  the  quantity  measured  by  means  of  nuclear  resonance  prof  i !  ir.g  is  the 
1  K 

"tal  0  concentration,  i.e.  /->  *  —  f  Q  *  n  +■  .  r>  *  \ 

P*,*-  ~  \Yl  *  Tr  +  /  »  I- 

" hi-  problem  of  solving  the  full  set  of  equations  will  not  he  addressed 
instead,  some  approximations  will  be  made  which  allow  great  simplif:- 
:>'U,  T!i<-  validity  of  the  approximation  may  of  course  be  quest  ieneu. 

■st  of  a.i  '  ,  it  can  be  seen  that  by  adding  equations  (D! )  through  (1H) 

dt 


he  approximation 


is  maiie,  then  can  be  found  by  finding  solutions  to  equations  (D4)  and 

■05),  which  thus  decouple  from  the  first  three. 


These  latter  equations  can  be  seen  to  revert  to  the  equations  generally 
:ken  to  describe  water  permeation,  which  assume  both  that  P* >y  Fa 

that  equilibrium  ? s  established  quickly  between  P  and  ^  .  Under  those 
on-lit. ions  adding  equations  (D4)  and  (.05)  results  in 

h  A  )  =  5  ’Vi  *  -k  (ft  -  £  )  . 


pnlying  the  inequality  as  well  as  the  equilibrium  condition  given  as  eq un¬ 
ion  (6)  in  the  test,  this  reduces  to 


l 

it 


y°r 


(l>&) 


0-1.  H.  W.  Ashcroft  and  N.  D.  Mermin,  Solid  State  Physics  (Holt  Rinehart,  1976) 
p  603 
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which  is  the  equation  that  both  Doremus  (2)  and  Warner  (3)  have. 


An  interesting  case  occurs  in  the  limit  that  T  is  so  large  that  the  mean 
free  path  between  water  dissociations  is  large  compared  to  the  SiCb  film 
thickness.  Under  those  conditions  A  approaches  its  saturation  value  while 

Pr  is  still  negligible.  The  buildup  of  pr  can  then  be  described  by  setting 
n  s  /°hi0  in  equation  (D5) .  That  equation  then  decouples  from  equa¬ 

tion  (D4)  as  well;  its  solution  is 


2jj±k£L  ($cir)  A 
P  oh  ($**' )  t 


independent  of  x. 


18 

Since  the  exchange  of  an  individual  0  between  the  network  and  dissolved 
molecular  water  depends  on  the  recombination  of  OH,  the  value  of  tracer  dif- 

fusivity  D*  at  any  time  is  proportional  to  the  value  of  at  that  time. 

Under  these  conditions,  the  apparent  value  of  D*t  at  any  time  is  related  to 
its  limiting  value  (i.e.  when  equilibrium  has  been  established  by  the  equation 


pet  reert- 


=  2> 


C  tt'  &  v  > 


This  equation  can  be  integrated  to  give 


I- 


Dl±  1 


1 t#sr\lr\ 


2-r  * 


o  .  A  *  ^ 

uWrt  K  = 


It  can  be  seen  that  the  quantity  [j3*0apparent  smaH  f°r  short  times;  it 

reaches  half  its  asymptotic  value  at  t  =  — .  Using  equation  (8)  in  the  text, 

this  time  is  approximate! v  100 t  at  600  C.  The  actual  value  of  x  at  any  tem¬ 
perature  can  be  determined  through  measurement  of  the  time  behavior  of 

r"n*n 

|_  ^apparent. 
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